RADC-TR-89-243 
Final  Technical  Report 

Nowmbar  1989 


AD-A217  035 


PRINCIPLES  OF  NONLINEAR  OPTICS 

Syracuse  University 
Partha  P.  Banerjee 


2  28 


This  report  has  been  reviewed  by  the  RADC  Public  Affairs  Division  (PA) 
and  is  releasable  to  the  National  Technical  Information  Services  (NTIS)  At 
NTIS  it  will  be  releasable  to  the  general  public.  Including  foreign  nations. 


RADC— TR— 89-243  has  been  reviewed  and  is  approved  for  publication. 


DOROTHY  JACKSON 


Project  Engineer 


>90HN  A.  GRANIERO 
Technical  Director 
Directorate  of  Communications 


FOR  THE  COMMANDER: 


BILLY  G.  OAKS 

Directorate  of  Plans  &  Programs 


If  your  address  has  changed  or  if  you  wish  to  be  removed  from  the  RADC 
mailing  list,  or  if  the  addressee  is  no  longer  employed  by  your 
organization,  please  notify  RADC  (DCLW  )  Grifflss  AFB  NY  13441-5700. 
This  will  assist  us  in  maintaining  a  current  mailing  list. 

Do  not  return  copies  of  this  report  unless  contractual  obligations  or 
notices  on  a  specific  document  require  that  it  be  returned. 


REPORT  DOCUMENTATION  PAGE 


form  Approved 
OMB  No.  07044188 


la.  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 

N/A _ 


2  b.  DECLASSIFICATION  /DOWNGRADING  SCHEDULE 

N/A 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 


6a.  NAME  OF  PERFORMING  ORGANIZATION 
Syracuse  University 


6c  ADDRESS  (City,  State,  and  ZIP  Code) 

Office  of  Sponsored  Research 
Skytop  Office  Building,  Skytop  Road 
Syracuse  NY  13210 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

Rome  Air  Development  Center 


ADORESS  (Gty,  State,  and  ZIP  Code) 


8b.  OFFICE  SYMBOL 
(If  applicable) 

DCLW 


tb.  RESTRICTIVE  MARKINGS 

N/A 


3.  DISTRIBUTION  /  AVAILABILITY  OF  REPORT 
Approved  for  public  release; 
distribution  unlimited. 


S.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

RADC-TR-89-243 


7a.  NAME  OF  MONITORING  ORGANIZATION 
Rome  Air  Development  Center  (DCLW) 


7b.  ADDRESS  (Orty,  State.  end/IPCode) 


Griffiss  AFB  NY  13441-5700 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

F30602-81-C-0169 


10.  SOURCE  OF  FUNDING  NUMBERS 


Griffiss  AFB  NY  13441-5700 


1 1 .  TITLE  (Include  Security  OasnticatJonJ 

PRINCIPLES  OF  NONLINEAR  OPTICS 


12.  PERSONAL  AUTHOR(S) 

Partha  P.  Banerjee 


13a.  TYPE  OF  REPORT  113b.  TIME  COVERED 

Final 


16.  SUPPLEMENTARY  NOTATION 


PROJECT 

TASK 

NO. 

NO 

2863 

92 

Mmm 1 1 


TO  Sep  88 


14.  DATE  OF  REPORT  (Yetr,  Month.  Oey)  115.  PAGE  COUNT 

November  1989  I  100 


COSATI  COOES 


GROUP  SUB-GROUP 


18.  SUBJECT  TERMS  (Continue  on  reverie  if  necessary  and  identify  by  block  number) 
Optical  Bistability 
Phase  Conjugation 
Harmonic  Generation 


19.  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

This  report  contains  a  summary  of  essential  principles  of  nonlinear  optics  such  as 
optical  bistability,  phase  conjugation,  and  harmonic  generation. 


20.  DISTRIBUTION  I  AVAILABILITY  OF  ABSTRACT 
C8  UNCLASSIFIEOflJNLIMITED  □  SAME  AS  RPT. 


22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Dorothy  Jackson 


00  Form  1473.  JUN  86 


121.  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


22b.  TELEPHONE  (Include  Are*  Code)  22c.  OFFICE  SYMBOL 

(315)  330-4092  RADC  (DCLW) 


Previous  editions  ere  obto/ete.  SECURITY  CLASSIFICATION  OF  THIS 

UNCLASSIFIED 


i 


CONTENTS 


Introduction  1 

2.  Origins  of  Nonlinearity  3 

3.  Harmonic/ subharmonic  Generation  11 

3.1  Mathematical  Formulation  for  One-dimensional  Propagation  13 

3.2  Harmonic  Generation  with  Quadratic  Nonlinearity  16 

3.3  Harmonic  Generation  with  Quadratic  and  Cubic  Nonlinearities  24 

4.  Self-refraction  28 

5.  Optical  Bistability  39 

5.1  Acoustooptic  Bistability  41 

5.2  The  Nonlinear  Fabri-Perot  Cavity  46 

5.3  The  Linear/Nonlinear  Interface  61 

6.  Phase  Conjugation  71 

6.1  Comparison  with  Holography  74 

6.2  Semiclassical  Analysis  77 

7.  The  Nonlinear  Schrodinger  Equation  and  Soliton  Propagation  81 


The  Nonlinear  Schrodinger  Equation  and  Soliton  Propagation  81 


1 .  Introduction 


This  presentation  on  nonlinear  optics  doe3  not  aim  to  be  complete  and 
detailed  in  the  sense  of  recalling  and  accounting  all  experiments  that  were 
performed  in  the  area,  and  carefully  presenting  rigorous  mathematical 
details.  Rather,  it  tries  to  present  the  essentials  of  the  physics  of 
nonlinear  optics  in  the  simplest  possible  way  and  unify  all  the  effects  under 
a  common  umbrella.  What  I  hope  to  achieve  is  a  simplified  treatment  of  most 
nonlinear  optical  effects  to  provide  the  inexperienced  reader  with  a 
springboard  to  jump  into  the  sea  of  existing  literature  which  is  often 
incommensurate  with  respect  to  notations  and  explanations. 

Without  further  ado,  let  me  briefly  describe  the  contents  of  this 
monograph.  Section  2  describes  the  origins  of  nonlinearity,  tracing  it  back 
to  its  manifestation  as  the  modification  of  the  linear  phase  velocity  of  a 
propagating  wave.  Relationships  are  established  between  the  nonlinearity 
parameters  we  choose  to  use  for  the  simplest  possible  description,  and 
existing  parameters  {seemingly  endless!)  in  literature.  MKS  unit3 
have  been  used  consistently.  Throughout  the  monograph,  z  is  assumed 
to  be  the  "nominal"  direction  of  propagation  and  a  "workhorse"  dependent 
variable  y  used  to  describe  the  physical  effect. 

Section  3  is  dedicated  to  describing  the  effects  of  harmonic/subharmonic 
generation.  One-dimensional  propagation  is  assumed  for  simplicity  and  the 
parametric  interaction  process  involved  is  described.  An  analysis  is 
advanced  for  harmonic/subharmonic  generation  in  a  quadratically  nonlinear 
medium  without  any  "dc"  polarization,  both  in  the  presence  and  absence  of 
dispersion.  It  i3  also  shown  that  similar  effects  may  be  observed  in  media 
with  quadratic  and  cubic  nonlinearities.  The  very  recent  observation  of 
harmonic  generation  in  cubically  nonlinear  dispersive  fibers  has  not  been 
treated  since  a  consolidated  theory  for  this  is  still  in  the  developmental 
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From  Section  4,  we  concentrate  solely  on  effects  caused  by  the  cubic 
nonlinearity.  In  this  section,  for  instance,  we  provide  a  simplified 
analysis  of  self-refraction  of  beams  which  occurs  due  to  an  induced 
refractive  index  profile  due  to  the  nature  of  the  propagating  beam.  Not 
mentioned  in  detail  is  the  phenomenon  of  "self-attraction"  between  adjacent 
beams,  although  a  figure  is  devoted  to  illustrating  the  effect.  More  recent 
topics  such  as  self -bending  of  an  asymmetric  beam  in  such  a  medium  is  left 
out  for  the  sake  of  simplicity. 

Another  very  interesting  nonlinear  optical  effect  is  bistability  and 
hyteresis,  which  is  addressed  in  Section  5.  To  illustrate  the  basic 
principles,  a  hybrid  acoustooptic  bistable  device  is  first  considered  and 
compared  to  the  electronic  Schmitt  trigger  to  highlight  the  necessity  of 
nonlinearity  and  positive  feedback  in  achieving  bistable  operation.  The 
nonlinear  Fabri-  Perot  cavity  is  next  analyzed  using  the  so-called  Maxwell- 
Bloch  equations  as  a  starting  point.  The  limiting  cases  of  purely  absorptive 
and  purely  dispersive  bistabilities  are  discussed  in  some  detail.  In  the 
latter  case,  the  onset  of  instability  is  addressed  by  considering  familiar 
feedback  criteria.  Lastly,  transmission  through  a  linear/nonlinear  interface 
is  analyzed  and  bistable  and  hysteretic  phenomena  predicted.  Bistable  optical 
elements  can  be  used  as  the  basic  building  blocks  of  an  optical  computer. 

The  presence  of  a  cubic  nonlinearity  can  be  sometimes  exploited  to 
achieve  real-time  holography  or  phase  conjugation.  This  is  discussed  in 
Section  6.  Principles  of  holography  are  recalled,  and  a  semiclassical 
treatment  of  phase  conjugation  (also  called  four-wave  mixing)  advanced  under 
the  simplifying  assumption  of  strong  "pumps." 

Finally,  in  Section  7,  the  role  of  a  cubic  nonlinearity  in  balancing 
dispersion  in  fibers  to  ensure  distortionless  pulse  propagation  is  analyzed 
through  the  nonlinear  Schrodinger  equation,  which  is  first  heuristically 
derived.  The  distortionless  pulses  arising  out  of  a  balance  between 
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nonlinearity  and  di3persion  are  commonly  called  solitons,  and  shew  promise 
of  future  use  in  fiber-optic  communication  systems. 

While  much  of  the  discussion  has  been  borrowed  and  "rebottled"  in  an 
unified  way,  seme  of  the  effects  reported  are  new  and  the  analysis  novel. 
Reference  may  be  made  to  the  idea  of  second  harmonic  generation  in  the 
presence  of  quadratic  and  cubic  nonlinearities,  the  comparison  of  the 
acoustooptic  bistable  device  with  an  electronic  Schmitt  trigger,  the 
treatment  of  purely  dispersive  bistability  in  a  Fabri-Perot  cavity,  and  the 
prediction  of  hysteresis  and  bistability  during  transmission  through  a  linear 
nonlinear  dispersive  interface. 

2 .  Origins  of  Nonlinearity 

In  an  optically  nonlinear  medium,  the  nonlinearity  can  be  attributed  to 
the  dependence  of  the  phase  velocity  on  the  amplitude  of  the  propagating 
wave,  unlike  the  linear  case  which  always  assumes  infinitesimal  wave 
amplitudes.  Thus,  while  a  PDE  for  a  wavefunction  y  of  the  form 

dy/5t  +  cQ3y/3z  -  0  (1) 


where  t  denotes  time  can  explain  unidirectional  propagation  along  z  in  a 

nondispersive  linear  medium,  the  phase  velocity  c  needs  to  be  modified  to  c 

o 

according  to 


P 


c 

P 


Co(1  +  &2V  +  ^3v2  +  -  •  '  > 


(2) 


to  account  for  nonlinear  propagation  [1] .  In  (2),  (3^  anc*  P3  are  referred  to 
as  the  quadratic  and  cubic  nonlinearity  coefficients  respectively  for  reasons 
that  will  shortly  become  clear.  The  wavefunction,  y  may  represent  a 
component  of  the  (real)  electric  field. 
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In  the  nonlinear  regime,  eq.  (1)  is  modified,  in  the  light  of  (2),  to 
3\y/3t  +  cq{1  +  P2v  +  P3y2)3y/3z  *>  0  (3a) 

or 

3y/3t  +  co3y/3z  +  (coP2/2)3y2/3z  +  ( c^P ^ / 3 )  3\y^/3z  =  0,  (3b) 

justifying  the  names  given  to  P2  and  P3  above.  We  remark,  chat  the 
nonlinearity,  e.g.,  P2,  is  responsible  for  shock  formation  during  fluid  flow 
due  to  steepening  of  the  (baseband)  wavefronts  as  shown  in  fig.  1.  This 
happens  since  the  parts  of  the  baseband  pulse  having  larger  amplitudes  travel 
faster  than  the  parts  with  smaller  amplitudes  for  P2>0,  leading  to  a  point  in 
time  where  the  right  edge  of  the  pulse  develops  infinite  steepness  or 
’shock.’  In  optics,  the  quadratic  nonlinearity  (and  cubic  nonlinearity,  too!) 
is  responsible  for  second  harmonic/subharmonic  generation,  and  thi3  will  be 
discussed  in  Section  3.  Special  effects  arising  from  the  cubic  nonlinearity 
will  be  mentioned  below  and  discussed  at  length  in  sections  4  -  7 . 

The  reader  may  realize  that  the  simple  model  described  above  cannot 
describe  wave  propagation  in  higher  dimensions,  necessitating  the  need  for  a 
higher  order  PDE.  This  can  be  readily  derived  from  (3)  by  differentiating 
w.r.t.  t  and  reusing  (3)  to  simplify.  Thi3  yields,  after  some  algebra, 

32y/3t2  -  co2aVaz2  ~  2co2  t  (P-/2)32v2/3z2  +  (p  /3)  d2y3/3z2]  (4a) 

s  P2dV/3t2  +  (2P3/3)3V/9t2  (4b) 

under  the  assumption  of  weak  nonlinearity  [2] .  This  assumption  helps  us 

neglect  second-order  terms  in  P  in  deriving  (4a),  and  in  approximating 

2/3 
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2  2  -2  2  2 

d  /9z  by  cq  9  /9t  on  the  RHS  of  (4a)  in  deriving  (4b) .  The  extension  to 

2  2 

higher  Tensions  may  be  effectively  done  by  replacing  9  / 9z  in  (4b)  by  the 


Laplacian 


V2: 


92y/9t2  -  cq2V%  -  3292V2/9t2  +  <2p3/3) 92y3/9t2 . 


(5) 


A3  will  become  clear  later  in  this  section,  eq.  (5)  with  P  t  0  becomes  the 

c.  t  6 

wave  equation  for  a  component  of  the  S  field  where  the  RHS  may  be  identified 
to  be  the  source  term  due  to  the  nonlinear  polarization  of  the  medium,  while 
the  quadratically  nonlinear  induced  polarization  is  responsible  for  second 
harmonic/ subharmonic  generation  as  stated  earlier,  the  cubically  nonlinear 
induced  polarization  gives  rise  to  an  amplitude-  dependent  refractive  index 
which  causes  self-refraction,  and  can  account  for  other  nonlinear  optical 
effects  such  as  bistability,  phase-conjugation,  soliton  propagation  etc., 
to  be  discussed  later. 

At  this  point,  it  is  instructive  to  strike  a  connection  between  the 
nonlinearity  parameters  P2,  P^  introduced  above  and  the  commonly  occurring 
parameters  in  nonlinear  optics  literature.  There,  the  description  of 
nonlinearitie3  is  given  in  terms  of  the  nonlinear  induced  polarization 
in  the  medium  or  an  amplitude  dependent  refractive  index.  In  the  former, 
the  induced  polarization  is  expressed  in  terms  of  a  Taylor  series 
expansion  of  the  electric  field  component  E,  according  to  [3] 


t  Z.jE.  +e  X.j.EjE,  +£  X.  ,E.E,  E,  +. 
oAij  3  o^ijk  j  k  o^ijkl  3  k  1 


(6a) 


A  e  X-'E-  +  2d.  ,,  E.E,  +e  X--.1E  •EuE1 
•  0*13  3  13k  3  k  o  i]kl  3  k  1 


(6B) 


A  e  x.  -ej  +  p. 

-  0*13  j  1 


NL 


(6c) 
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where  the  Einstein  convention  for  summation  over  repeated  indices  ha3  been 
employed.  In  (6),  is  the  linear  susceptibility,  while  Xj_jk  and  X^j^ 

refer  to  the  second  and  third-order  nonlinear  optical  susceptibilities 
respectively .  The  reader,  not  quite  familiar  with  tensor  notation,  ha3 
scope  for  being  baffled  for  sometime  at  first,  but  may  take  heart  from  the 
fact  that  (6)  is  merely  a  shorthand  for  writing  lengthy  expressions.  For 
example,  without  the  last  two  terms  of  the  RHS  of  (6),  P  ,  the  x-component 
of  the  polarization  may  be  expressed  as 

P1  ■  eo<XllEl  *  *12E2  +  *13E3'  <7» 

where  E^'E2,E3  stand  for  the  x,y  and  z  components  of  the  electric  field  in  a 
medium  where,  in  general,  xil  *  X12  *  ^13'  the  medium  anistropic. 
Readers  may  be  reassured  that  the  new  representation  is  not  usually 
necessary  to  explain  the  nonlinear  effects  mentioned  earlier,  but  should  be 
familiar  with  it  in  case  of  cross-referencing  with  existing  literature. 

The  wave  equation  for  the  vector  electric  field  E  may  be  readily 
derived  from  Maxwell's  equations  by  writing 

5  -  e  E  +  P  (8) 

o 

where  the  components  of  P  are  specified  in  (6) .  This  yields 

d2E/dt2  -  c  2V2i  -  -  (H  /e)d2pNL/dt2  (9) 

o  o 

in  a  region  free  of  sources  and  currents  and  for  a  medium  which  is  isotropic 
in  the  linear  regime  <X. .  “  XT) •  In  (9), 

li  L 


e  -  eQ(l  +  xL> 


(10a) 
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and 


c 

o 


<l/Hoe) 


1/2 


(10b) 


It  is  true  that  through  eq.  (9) ,  one  can  obtain  a  set  of  coupled  PDEs, 

in  general,  that  describes  the  behavior  of  a  certain  component  of  the 

electric  field.  However,  as  will  be  shown  later,  similar  equations  may  also 

be  derived  by  starting  from  (5)  and  assuming  no  particular  anisotropicity  in 

the  medium.  To  explain  the  physics  of  what  is  happening,  we  shall  therefore 

use  eqns .  (3)  or  (5)  as  our  starting  point. 

To  obtain  typical  values  for  for  instance,  note  that  a  comparison 

of  (4b)  and  (9),  with  (6),  shows  that  is  of  the  order  of  -2  (Jl^/e ) d_^  = 

-2u.  e  (1  +  X  )  d.  .  Listed  below  are  some  typical  values  for  d\  along 
*0  o  13k  13k 

with  the  corresponding  wavelengths  at  which  they  are  measured  [3] : 


crystal 

— 

X  ()lm) 

1  -22  3  2 

x  10  m  /V-s  ) 

GaAs 

10.6 

d123 

:  107 

InSb 

28 

dl23 

:  462 

CdTe 

28 

d123 

:  4  8 

CdS 

10.6 

d333 

:  35 

j 

d311 

:  -  21 

1 

1.06 

d333 

:  80 

LiNb03 

1.06 

d333 

:  -  27 

d222 

:  4 

quartz 

1.064 

dlll 

:  0.4 

KDP 

0.6328 

d213 

:  0.57 

1.06 

d312 

:  0.50 

Table  1.  Quadratic  nonlinearity  coefficients  d^^  ^or  some  crystals. 
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We  caution  readers  that  the  values  of  d, ,,  in  Table  1  are  only  valid 

ijk 


for  the  case  of  second  harmonic  generation.  For  the  general  three-frequency 

o  ■  (0  +03^),  the  values  of  d. 

O  i  A-  1  J  JC 


interaction  case  (ox,  *  (i^+ox,),  the  values  of  d^^  may  be  different  than 


those  listed. 

Now  to  prescribe  typical  values  for  the  cubic  nonlinearity 

coefficient,  we  will  first  state  its  relationship  with  the  amplitude- 
dependent  refractive  index  coefficient  n^  defined  as  [3] 


n  =  n  +  (1/2) n_y  y 
o  2  e  6 


(id 


where  represents  the  envelope  of  y:  =  Re  exp  j  (co^t-k^z) } ,  and  where 

n^  denotes  the  linear  refractive  index.  This  can  be  readily  achieved  by 
starting  from  (2)  with  (3 2  =  0  and  invoking  the  weak  nonlinearity  assumption. 
The  result  i3 


B_  =  -2n„/n  . 
k3  2  0 


(12) 


2  2 

The  units  of  n^  and  (3^  are  m  /v  in  MKS  since  has  units  of  V/m.  We 
specifically  point  thi3  out  for  the  benefit  of  readers  who,  like  us,  will 
no  doubt  be  frustrated  at  different  unit  systems  that  are  used  in 
literature,  as  well  a3  different  definitions  of  the  nonlinear  polarization 
and  the  nonlinear  refractive  index  coefficient.  For  instance,  an  equivalent 
representation  of  the  nonlinear  refractive  index  that  appears  in  literature 
is  [4] 


n  =  n  +  VI, 
o 


(13) 


where  I  is  the  optical  intensity  or  irradiance.  y  has  the  unit  of  m  / W  and  is 
related  to  n^  as 
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7  -  n„/c  e  n  .  (14) 

'  2  o  o  o 

Finally,  is  related  to  a  nonlinear  susceptibility  viz.  33 

n2  -  3Zuil/no,  (15) 

for  the  case  of  an  electric  field  linearly  polarized  in,  say,  the  x- 

NL 

direction.  Eq.  (15)  may  be  derived  by  writing  E  and  P  in  (9)  m  terms  o: 

NL  NL  2  * 

their  slowly  varying  envelopes  E  and  P  ,  setting  P  -  3E  Y....E  E  (3ee 
J  J  ’  e  e  e  0*1111  e  e 

(6b))  and  finding  the  effective  propagation  constant.  The  factor  3  above  is 
a  degeneracy  factor  which  arises  in  the  interaction  process  (0^  =  (O^+U^-co^ . 
Values  of  ^or  different  materials  are  listed  in  Table  2  [5]  . 


— 

Material 

X(fim) 

n 

0 

,4 Jt  ,  -21  2  ,„2> 

Xllll <_ 9  X  10  m  /V  ) 

! 

GaAs 

10.6 

3.3 

120 

InSb 

5.3 

4 

-6x10 10 

10.6 

4 

2xl06 

CdTe 

1.06 

3 

2.5xl05 

CdS 

2.42 

130 

G© 

10.6 

4 

io3 

Si 

10.6 

3.4 

60 

H2° 

0.694 

1.33 

°-7 

Acetone 

0.694 

1.35 

1.8 

Penzene 

1.064 

1.5 

2.4 

Table  2.  Cubic  nonlinearity  coefficients  %  for  some  materials. 
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3 .  Harmonic/Subharmonic  Genera1- ion 

An  important  consequence  of  nonlinearity  i3  frequency  multiplication. 

To  see  what  this  means  consider,  frrst,  one  dimensional  CW  propagation  in  a 

quadratically  nonlinear  (p^  0,  (3^  =  0)  nondispersive  medium.  If  a  CW  wave 

of  angular  frequency  co  and  propagation  constant  k  is  incident  on  such  a 

medium,  it  can  generate  harmonics  that  travel  along  with  the  same  velocity. 

The  energy  and  momentum  conservation  law3  as  well  as  the  dispersion 

relationship  (variation  of  co  with  k)  have  to  be  satisfied  whenever  wave 

interactions  occur  [6].  Limiting  ourselves  to  only  two  frequencies,  namely 

the  fundamental  (co  ,  k  )  and  the  second  harmonic  (2co  ,  2k  )  ,  the  following 

o  o  o  o 

interaction  occurs: 


CO  +  CO  ->  2co  ;  2co  -  CO  -»  CO 
o  o  o  o  o  o 

(16) 

k  +  k  -4  2k  ;  2k  -  k  -»  k 
o  o  o  o  o  o 

with  oWk^  =  cq/  the  phase  velocity  of  the  two  propagating  frequencies.  The 

relation  (15)  depicts  explicitly,  for  example,  the  paranetric  mixing  of  two 

waves  at  (co^k^)  through  a  quadratic  nonlinearity  to  yield  a  third  wave  at 

(2co  ,  2k  ) .  This  is  the  basis  of  second  harmonic  generation  usually  observed 
o  o 

in  nonlinear  electrooptic  crystals  [3],  [4],  [7], 

The  next  question  to  ask  is:  Could  a  cubic  nonlinearity  also  play  a 
role  in  second  harmonic  generation?  At  first  thought,  this  seems  doubtful, 
since  a  cubic  nonlinearity  can  only  facilitate  interactions  of  tne  type  d i)q  + 
C0q  +  C0q  -»  3(0^  leading  to  third  harmonic  generation.  A  closer  examination  has 
revealed,  however,  that  the  answer  is  yes,  but  only  if  either  the  quadratic 
nonlinearity  is  also  present,  or  if  a  "dc"  source  of  energy  is  either 
externally  applied  or  internally  generated  [8] .  The  latter  case  is  in 
support  of  recent  observations  of  second  harmonic  generation  in  optical 
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fibers  which  can  only  have  odd-order  nonlinearities .  We  will  not  treat  this 
case  since,  at  this  time,  no  undisputed  exp' anation  is  available.  In  the 
former  case,  the  parametric  mixing  that  occurs  can  be  represented  as: 

2(0  +  (0  -  CO  ->  2(0  ,  2(0  -  2(0  +  2'J  ^  2(i)  ; 

ooo  oo  O  0  o 

(0  +2(0  -2(0  -4(0,(0  -  (0  +  (0  ->  a>  ,  (17) 

o  o  oooooo 

with  similar  relations  for  the  k's.  The  parametric  process  described  above 

may  be  visualized  as  contributic  a  to  the  growth  of  the  second  harmonic 

amplitude  in  the  presence  of  tta  fundamental  (and  vice-ver3a) . 

If,  now,  the  medium  under  consideration  is  dispersive,  the  propagation 

constants  for  the  frequencies  (0q  and  2o)q  are  no  longer  related  by  a  simple 
A  A 

factor  (viz.,  k ^  =  k(2a)Q)  =  2k^  =  k((0Q)).  For  one-dimen3ional  propagation, 
this  gives  rise  to  a  spatial  beating  for  the  fundamental  and  harmonic 
(subharmonic)  amplitudes,  which  in  the  general  theory  of  nonlinear  waves,  i3 
called  the  "Fermi-Past* -Ulam”  recurrence  phenomenon  [9].  The  recurrence 
period,  which  in  the  linear  case  for  two  propagating  frequencies  is  given  by 
2n/ (k2  —  2k^) ,  is  mor if ied  by  the  presence  of  nonlinearities,  to  be  explained 
in  detail  below. 

In  many  physical  situations,  however,  wave  propagation  is  essentially 
higher-dimensional  unless  special  care  is  taken  to  ensure  one-dimensional 
propagation.  For  the  higher-dimensional  dispersive  case,  the  two 
frequencies  are  no  longer  constrained  to  propagate  in  the  3ame  direction, 
and  hence  interact  according  to  the  "resonant-triad"  relations  [6] 


(0i  +  (02  +  co3 


kl  +  k2  -  k3 


(18a) 
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with 

(0  =  CO(  |k^ ! )  .  (18b) 

The  resonance  triad  vector  diagram  is  shown  in  fig.  2  for  the  case  co  =  to 

(0  and  <D_  -  2(0  .  Since  the  waves  are  "free"  to  travel  in  higher  dimensions 
o  3  o 

(predictable  from  the  vector  diagram  and  the  dispersion  relationship),  there 
is  no  longer  any  spatial  beating  due  to  dispersion,  and  the  evolution  of  the 
second  harmonic/subharmonic  is  similar  to  the  nondispersive  case,  albeit  in 
their  respective  directions  [6],  [10],  [11].  The  rigorous  mathematical 

formulation  of  this  case  is  beyond  the  scope  of  this  discussion. 

3 . 1  Mathematical  Formulation  for  One-dimensional  Propagation 
We  first  assume  y(z,t)  to  be  representable  as 

2 

y(z,t)  =  -  y~  ¥  (z)exp  [  jn  (CO  t  -  k  z)  ]  (19) 

2  i — _  n  o  o 

n=-2 

★ 

with  ¥  =  0  (no  dc  component) ;  and  with  V  =  'F  to  ensure  that  the 
o  -n  n 

wavefunction  is  real.  Then  substituting  (19)  into  (3)  and  gathering  terms 

around  the  respective  frequency  components,  viz.,  (0  and  2co  ,  we  obtain 

o  o 

the  following  3et  of  coupled  spectral  evolution  equations  [8] : 

2 

dV  /dz  -  (jk  p./2)  r  mV  V 

n  o  2  z_,  m  n-m 


+  'VF”  i  i 

m»-2  l*-; 


m 


(20) 


In  deriving  (20),  we  have  assumed  (o  /k  -  c  and  that  *F  (z)  is  a  slowly 

o  o  o  n  J 
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u>  1='^2=  1,0  3'/2 


Iks  I 

2IM 


Fig.  2.  Resonant  triad  wavevector  diagram  for  three  interacting  waves 

1,  2  and  3,  two  of  which  have  the  same  temporal  frequency.  The 
wavevector  diagram  follows  from  the  conservation  of  momentum 
during  wave  interactions  and  the  dispersion  characteristics  of 
the  system.  The  interaction  angle  9  depends  on  the  amount  of 
dispersion. 
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varying  function  of  z  in  the  sense  that  Id'P  /dzl  «  I  nk  I  .  Note  that  the 
first  and  second  terms  on  the  RHS  of  (20)  respectively  denote  the 
contributions  from  the  quadratic  and  cubic  nonlinearities. 

To  account  for  the  effect  of  dispersion,  we  must  3lightly  modify  (20) . 
To  do  this,  we  may  realize  that  in  the  presence  of  dispersion  alone,  the 
propagation  constant  for  the  'n'-th  har  monic  can  be  written  as  k  .  Thus, 
during  propagation,  the  spatial  behavior  of  4^  will  be  [12] 

4*  =  V  (0)  exp [- j (k  -nk  )z],  (21) 

n  n  no 


so  that  in  (19)  , 


^(z,t) 


2 

d/2)  y 

n=*-Z' 


4*  (0)  exp  [  j  (n<0  t-k  z)], 
n  on 


as  expected.  Now  (21)  is  analogous  to 


(22) 


d'F  /dz  =*  -j  (k  -nk  )  4*  .  (23) 

n  non 

To  incorporate  the  effect  of  dispersion  in  our  nonlinear  system  (20),  we  may 
heuri3tically  add  a  term  as  on  the  RHS  of  (23)  to  the  RHS  of  (20)  if  the 
nonlinear  and  dispersive  effects  are  small: 

2 

d4*  /dz  -  (jk  (3_/2)  y  m4/  4/ 
n  J  or 2  L-j-  ra  n-m 

m=-2 


+  ( jJcQ03/4) 


4/ 

m-£ 


4> 

n-m 


-j  (k  -nk  )  4*  . 
n  on 


(24) 
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Eq.  (24)  will  serve  as  our  model  equation  describing  the  evolution  of  the 
fundamental  and  the  second  harmonic  in  a  medium  with  quadratic  and  cubic 
nonlinearities,  and  dispersion.  In  what  follows,  we  will  consider  two 
special  cases:  (a)  harmonic  generation  with  quadratic  nonlinearity  without 
and  with  dispersion,  and  (b)  harmonic  generation  with  quadratic  and  cubic 
nonlinearities  and  no  dispersion.  Comments  on  subharmonic  generation  will 
also  be  made . 

3 . 2  Harmonic  Generation  with  Quadratic  Nonlinearity 

In  this  case,  with  *  0/  P3  =  0,  eq.  (24)  gives  the  following  set  of 
equations  for  n  -  1  and  2: 

dV./dz  -  j(P-k  /2)'P.'P.  -  j(k  -k  )¥.,  (25a) 

l  2  o  —  12  1  o  •  1 

d*P2/dz  -  j(P2kQ/2)‘i'12-  j(k2-2ko)«?2,  (25b) 

where  VF_ ,  'P_  denote  the  complex  spectral  amplitudes  at  co  and  2(0  .  Note 
12  OO 

that  we  have  retained  the  symbols  k^  and  k2  to  3tand  for  the  propagation 

constants  at  the  frequencies  (0q  and  2(0q  respectively  at  this  point.  In 

principle,  k.  could  be  set  equal  to  k  without  loss  of  generality;  however, 
1  o 

we  choose  to  retain  the  distinction  to  bring  out  the  effect  of  dispersion 
more  succinctly.  We  will  now  resolve  the  complex  spectral  amplitudes  into 
their  respective  (real)  amplitudes  and  phase  by  writing 

*P  -  a  exp(-j<b  ]  -  'P  (26) 

n  n  n  -n 

in  (25).  This  yields  the  following  set  of  four  equations: 


da^/dz  -  (P2ko/2)  a1a2  sin  (<t>2~24> ^ )  ,  (27a) 

da2/az  -  -  (02ko/2)  sin  (<l>2-24> x )  ,  (27b) 

d^/dz  =  -(P2kQ/2)a2  cos  id  ^  — 2<}>  ^ )  +  (k^-k^),  (28a) 

dty  /dz  =  -(P-k  /2)  (a2/a  )  cos (0  -2$  )  +  (k  -2k  ).  (28b) 

Z  Z  O  Z  Z  1  Z  O 


Now,  (28a, b)  can  be  combined  by  defining 


0  =  <t»2  — 2<D  x  (29) 


to  give 


d0/dz  *•  Ak  +  ( P2 k^ / 2 )  (2a2~a2/a2)  cos  0,  (30) 


where 


Ak 


A 


(31) 


i3  proportional  to  what  is  often  called  the  "phase  velocity  mismatch"  due  to 
dispersion.  Remark  that  from  (27a, b),  it  readily  follows  that 


2  2 
ai  +  a2 


constant 


a2  (0) 


a2(0) 


E, 


(32) 


depicting  conservation  of  energy.  Finally,  normalization  using  the 


definitions 
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-1/2 

(0,k  /2)EX  z, 
z  o 


As  -  2Ak/  (jJlc^E172) , 


,~l/2 
a  /E  , 


reduces  the  system  (27),  (32)  to 


du/d£  -  uv  sin  0, 


2 

dv/d£  «  -u  sin  0, 


d0/di;  ■  As  +  (2v-u2/v)  sin  0, 


(33) 


(34a) 


(34b) 


(35) 


where  we  have  used  the  definition  of  8  as  in  (29) .  Observe,  also,  that 
using  (34a, b) , 

d(fcn  (u2v)  ]  /d^  -  (2/u)du/dJj  +  (l/v)dv/dij 

=  (2v-u2/v)sin  8;  (36) 

hence,  (35)  may  be  recast  in  the  form 

d0/d£  -  As  +  cot  0  d[in (u2v) ] /d£ .  (37) 


Eqs.  (34),  (37)  therefore  describe  the  evolution  of  the  real  amplitudes  and 
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the  relative  phase  difference  during  propagation  in  a  quadratically 
nonlinear  dispersive  medium  [7],  [13]. 

Consider,  first,  the  case  when  As  -  0  (perfect  phase-matching) .  Then,  from 
(37) ,  one  integration  yields 


2 

u  v 


cos 


constant 


S  r  = 


u  (0 ) v (0) cos0  (0) 


(38) 


Now,  using  (34b),  (38)  and  realizing  that 


2 

u 


2 

+  v 


1, 


(39) 


we  can  write 


?  2  22221/2 
d(v  )/d^  =  *  2v(l-v  )  (1  -  TZ/v  U-vV] 


or 


\  =  ±(1/2)  f  ^  [v2(l-v2)2-  r2]~1/2  d(v2 

v  (0) 


(40) 


The  general  solution  to  v(Jj)  (and  hence  to  u(£)  and  0(£)),  as  given  by 
(40),  is  rather  complicated  if  v(0)  *  0.  Suffice  here  to  state  that  the 
solution  depends  on  the  roots  of  the  equation 


2„  2^2  r2  n 

v  (i-v  )  -  r  -  o 


(41) 


2  2  2 

which  may,  for  now,  be  written  as  v  ,  v^  and  v  ,  with 

a  b  c 


2  2  2 
V  >  V,  >  V  . 

c  —  b  —  a 


(42) 


The  solution  is  expressible  in  terms  of  Jacobian  elliptic  functions  like 
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'on'  and  'an*  which  look  similar  to  the  'cos'  and  'sin'  functions,  and  are 
periodic  in  nature,  with  a  period  IX  given  as  [14] 


2  2  2 
[v  (1-v  ) 


r2  -1/2  2 

r  ]  d (v  ) . 


(43) 


A  case  of  physical  interest  is  where  v(Q)  -  0  (u(0)  *0),  i.e.,  when  we 

want  to  generate  the  second  harmonic  starting  from  the  fundamental.  Then  T  - 

2  2  2 

0  (from  (38))  and  v^  **  0,  *  v^  =  1  (from  42)) .  The  period  (using  (43)) 

becomes  infinity.  (40)  may  be  readily  integrated  to  give 


v(£)  =  tanh  § 


(44) 


so  that  from  (39),  a  possible  solution  to  u  is 

u<5)  “  sech  J; .  (45) 

Now,  (44),  (45)  suggest  that  at  =  0,  8(0)  =  ~n/2  (using  (34))  =  8(^)  (using 

(35) )  since  As  ■  0  (perfect  phase  matching) .  A  plot  of  u  and  v  is  shown  in 
fig.  3a.  Note  that  the  fundamental  amplitude  asymptotically  goes  to  zero; 
while  the  second  harmonic  amplitude  asymptotically  tends  to  1,  although  in  a 
physical  system  with  attenuation,  the  second  harmonic  amplitude  eventually 
decays  to  zero. 

We  remark,  in  passing,  that  eq3 .  (34),  (35)  may  also  be  used  to  explain 

subharmonic  generation.  For  this,  we  have  to  take  v(0)  *  0,  and  assume  a 
small  amount  of  the  subharmonic  (u)  exists  in  the  system  due  to  ambient 
noise.  If  we  assume  this  initial  amount  of  the  subharmonic  (-ev(O)),  it  is 
easy  to  check  that  for  9(0)  -  it/ 2,  the  solution  to  (34)  may  be  expressed  as 

u(£)  -  sech  (£  -  4q) ,  (46a) 


(46b) 


v($)  «  -tanh  (^  -  §  )  , 


with 


-  tanh  1  v (0 )  r  (1/v (0 ) ) £n ( 1/e ) 

and  where 

V2  (0)  (1+e2)  =  i, 

to  ensure  conservation  of  energy.  A  plot  of  u  and  v  for  this  case  i 
in  fig.  3b. 

Consider,  next,  the  case  when  As  #  0  (imperfect  phase-matching) 

2 

this  case,  multiplying  (35)  with  u  v  sin  0  and  using  (34),  we  get  [7 
(d/d£) (u2v  cos  0]  +  (As/2)  d(v2)/d£  -  0. 

Integration  gives 

2  2  A 

u  v  cos  0  +■  (As/2)v  =  constant  =  T. 

As 

“  T  +  (As/2)  v2  (0)  . 

Eq.  (40)  is  now  generalized  as 

v2  (£) 

%  -  ±  (1/2)  [  2  [V2 (1-V2) 2  -  {r  -  (As/2) (v2-v2(0))}2]'1/2d(v: 

Jv  (0) 

Everything  previously  said  about  (40)  and  its  solutions  for  v(0)  * 


(47) 

(48) 

s  shown 

.  In 
],  [13] 

(49) 

(50) 

V 

(51) 

0  carries 
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over  to  the  solution  of  (51),  and  the  period  in  this  case  is  given  by 
2 

V  „  ^ 

-  j  2b  [v2(l-vV  -  {r  -  (As/2)  (v"-v2(0))}2]"1/2  d(v2),  (52) 

Va 


where  v  ,  v,  ,  v  (with  v  >  v,  >  v  )  denotes  tne  roots  of  the  quantity  in 

C  JO  3  0  0  3 

square  brackets  in  the  above  expression. 

As  before,  a  case  of  physical  interest  is  where  v(0)  =  0  (u(0)*0).  We 
will  find  the  period  of  oscillation  for  this  case.  To  do  this,  we  recast 
(52)  into  the  form 


(2/v  ) 
c 


71/ 2 


2  2 

da  [  1-  (v.  /v  )  sin  a] 
o  c 


■1/2 


(53) 


by  employing  the  substitution 


.  2 

sin  a 


,  2  2  .  2  2X 
(V  -V  )  /  (V  -V  ) 

a  b  a 


(54) 


and  noting  that  for  v(0)  =0,  T  *  0  and  v  =0.  The  other  two  roots  are 

d. 

given  by  the  solution  to  the  algebraic  equation 


2  2  2  2  2 
v  (1-v  )  -  {(As/2)v  }  =  0 


(55) 


as 


2  2  1/2  2  2 

vq  -  (As/ 4)  +  { 1+  (As/ 4)  }  ]  =>  1/Vfc. 


(56) 


For  the  severely  mismatched  case  (As>>1),  it  can  be  shown,  after  simple 


algebra,  that 


v2  =  (2/Aa)2  [1-8/ (As)2] 
b 
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(57a) 


v2  ~  (As/2) 2  [l+2/[ As] 2], 
c 


(57b) 


[The  period  in  z  may  be  readily  found  using  the  tran3f ormations  in 
(33)].  The  growth  of  the  second  harmonic  amplitude  for  different  values 
of  phase  mismatch  is  shown  in  fig.  4.  In  the  theory  of  nonlinear  waves, 
this  periodic  behavior  is  referred  to  as  the  Fermi -Pasta-Ulam  (FPU) 
recurrence  [9]  . 


3 . 3  Harmonic  Generation  with  Quadratic  and  Cubic  Nonlinearities 

We  will  discuss,  in  this  section,  the  problem  of  harmonic  generation 
with  quadratic  and  cubic  nonlinearities,  and  no  dispersion.  The  case  with 
dispersion  follows  in  a  similar  way,  and  will  not  be  discu33ed  here,  since 
our  intention  is  to  simply  point  out  the  criteria  that  dictate  harmonic 
generation  in  a  quadratically  and  cubically  nonlinear  environment. 

For  the  analysis,  we  return  to  eq.(24),  and  neglect  dispersion  to  write 
down  the  following  set  of  equations  for  the  spectral  amplitudes  and  4^ 

[8]  : 

d^/dz  -  j(hQP2/2)  4'_14'2  +  j(k  0  /4)  ('P_1'l£  +  2'VJ?2'V_2)  (59a) 

d^/dz  -  j(koP2/2)^  +  j(kQP3/2)  tf2*_2  +  V'-l'*2>*  (59b) 


Now,  using  the  same  procedures  as  that  followed  in  Section  3.2,  we  can 


o 


Fig.  4.  Growth  of  the  second  harmonic  amplitude  with  propagation  in  the 
phase-mismatched  case.  For  comparison,  the  phase-matched  case 
is  included  in  the  figure.  The  approximate  expression  for  the 
period  is  given  in  eq.  (58). 
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derive  the  coupled  equations  between  the  (normalized)  real  amplitudes  and 
the  relative  phase  difference.  Relations  to  use  are  (26),  (29),  (32)  and 
(33)  without  the  phase  mismatch  term  As,  and  the  definition 

Q  -  <P3/232)  E1/2.  (60) 

The  quantity  Q  may  be  thought  of  as  a  parameter  describing  the  ratio  of  the 
cubic  and  quadratic  nonlinearity  coefficients.  The  set  of  coupled  equations 
for  this  case  eventually  becomes 

du/d!j  •  uv  sin  0,  (61a) 

2 

dv/dg  »  -u  sin  0  (61b) 


and 

d0/df;  ■  2Q(v2-u2)  +  (2v-u2/v)cos  0.  (62) 

As  expected,  eqs .  (61),  (62)  for  Q  =  0  (no  cubic  nonlinearity)  become 

identical  to  eqs.  (34),  (35)  for  As  *  0  (no  dispersion). 

To  solve  the  above  system,  we  proceed  analogous  to  the  method  in 

2 

section  (3.2)  for  As  *  0 .  For  instance,  multiplying  (62)  by  u  v  sin  0 
and  using  (61)  to  simplify,  we  get,  after  some  algebra, 

2  _  2  2  ^  2  2  2 

u  v  cos  0  -  Qu  v  ■  constant  -  T  -  u  (0)v(0)  cos  0(0)  +  Qu  (0)v  (0) . 

o 

(63) 


Finally,  using  (61b),  (39)  and  (63),  and  integrating,  we  obtain 
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4 


±(1/2) 


fVVa-vV  -  {r  -sa-v2)v2)2j-l/2  div2) . 
Jv  (0) 


(64) 


The  quantity  in  square  brackets,  unlike  the  similar  expressions  in  (40)  and 

2 

(51)  have  four  roots  since  it  is  a  quartic  polynomial  in  v  .  However,  it  is 
possible  to  reduce  (64)  to  a  problem  where  the  denominator  is  a  cubic 
polynomial,  through  a  simple  change  of  variable,  but  that  is  out  of  the  scope 
of  this  discussion.  Suffice  to  say  that  once  this  is  achieved,  everything 
previously  said  about  (40)  and  its  solutions  may  apply.  In  what  follows,  we 
will  consider  the  case  of  physical  interest  where  v(0)  =  0  (u(0)*0);  i.e., 
when  we  want  to  generate  the  second  harmonic  starting  from  the  fundamental. 
Then  rQ  =  0  and  (64,  becomes 

2 

-  fV  2  2  2  -22  -1/2  2 

5  -  ±<1/2Q)  [v  (1-v  )  (1/Q  -v  ) ]  d (v  ) .  (65) 

J0 


-2 

For  Q  >  1,  direct  integration. 


after  some  algebra,  gives 


v2(^)  -  1/[Q2  +  (Q2— 1 ) cot2  { (Q2-l) ] . 


(66) 


The  fundamental  and  second  harmonic  amplitudes  vary  periodically  with 
Jj,  with  the  period  given  by 


n/<52-l>1/2 


(67) 


This  indicates  recurrent  behavior,  similar  to  the  case  of  harmonic 
generation  with  quadratic  nonlinearity  and  in  the  presence  of  dispersion. 
All  nonlinear  materials  have  quadratic  as  well  as  cubic  nonlinearities,  and 
the  simple  analysis  above  tells  us  that  periodic  exchange  of  energy  will 
occur  if  the  condition 
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P3/2p2u(0)  >  1  (68) 

is  satisfied.  The  variation  of  the  fundamental  and  second  harmonic 

amplitudes  for  this  case  is  shown  in  fig.  5. 

~2 

On  the  contrary,  for  Q  <  1,  direct  integration  of  (65)  gives 

2  —2  —2  2  —21/2 
V  ($)  =  [Q  +  (1-Q  ) coth  { (1-Q  )  ] .  (69) 

This  suggests  an  asymptotic  increase  in  the  second  harmonic  amplitude  with 
graduate  decay  of  the  fundamental  as  shown  in  fig.  6,  and  is  similar  to  the 
case  of  harmonic  generation  in  the  presence  of  quadratic  nonlinearity  and  no 
dispersion . 

4 .  Self-refraction 

As  one  of  the  first  effects  of  a  cubic  nonlinearity  on  wave  propagation, 
we  will  first  consider  self-refraction  of  a  beam.  To  get  a  physical  picture, 
consider  a  beam  with  a  Gaussian-like  profile  as  shown  in  fig.  7.  If  the 
medium  i3  cubically  nonlinear,  we  can  describe  the  nonlinearity  in  terms  of 
n^,  the  nonlinear  coefficient  of  the  refractive  index,  which  is  related  to  P^ 
as  given  in  (12)  .  If  n^  is  greater  than  zero,  regions  along  propagation  with 
greater  amplitude  will  possess  a  greater  refractive  index  than  a  region  with 
lower  amplitudes.  The  result  is  a  refractive  index  profile  much  like  graded 
index  optical  fibers.  If  we  trace  'rays',  these  would  therefore  appear  to 
bend  towards  the  axis  of  propagation,  indicating  a  reduction  of  the  beam 
waist  size,  and  hence  an  increase  in  the  on-axis  amplitude.  This  simple 
picture  would  however  suggest  that  the  on-axis  amplitude  should  tend  to 
infinity;  however  this  does  not  occur  since  diffraction  puts  a  limit  to  the 
minimum  waist  size.  This  is  shown  in  fig.  7.  Heuristically  speaking,  this 
makes  sense  since  the  amount  of  diffraction  (as  predictable  from  the  angle  of 
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Fig.  5.  Variation  of  the  fundamental  and  the  second  harmonic  with 

propagation  in  the  presence  of  quadratic  and  cubic  nonlinearities 
and  no  dispersion.  The  parameter  Q  =  3.5. 


Variation  of  the  fundamental  and  the  second  harmonic  with 
propagation  in  a  medium  with  quadratic  and  cubic  nonlinearities 
and  no  dispersion.  The  parameter  Q  <  1. 
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(a) 


(b) 

Fig.  7.  Self-focusing  of  a  beam  in  a  medium  with  a  negative  cubic 

nonlinearityCB^ <  0) .  (a)  ray  trajectory:  rays  tend  to  converge 
toward  the  direction  of  maximum  refractive  index;  (b)  periodic 
focusing.  This  occurs  since  the  beam  waist  can  decrease  only 
till  the  point  nonlinear  effects  dominate  over  diffraction. 
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diffraction)  depends  on  the  ratio  of  the  wavelength  to  the  waist  size.  For 
an  arbitrary  initial  beam  profile,  therefore,  we  would  expect  initial 
reduction  of  beam  waist  size  before  diffraction  effects  start  to  dominate  and 
spread  the  beam  again,  resulting  in  periodic  focussing.  It  turns  out  that 
while  this  is  mostly  true  for  a  beam  in  a  two-dimensional  geometry,  it  may 
not  be  true  for  the  three-dimensional  case.  It  is  also  possible  to  find  the 
right  beam  profile  for  which  the  beam-narrowing  effect  of  nonlinearity 
exactly  balances  the  beam-spreading  effect  of  diffraction  [4],  [15],  [16], 
[17];  this  will  be  derived  shortly. 

For  the  opposite  kind  of  nonlinearity,  it  is  easy  to  argue  that  the 
beam  will  spread  more  than  it  does  for  the  linear  diffraction-limited  case 
(see  fig.  8)  . 

In  what  follows,  we  will  first  derive  an  evolution  equation  for  an 
arbitrary  beam  profile  in  a  cubically  nonlinear  medium  and  derive  an 
analytic  expression  for  the  diffraction-free  beam  in  the  medium.  We  write  y 
as 


V*  (1/2) Ye (x, y, z)  exp  j (O^t-k^z)  +  c.c.  (70) 


and  substitute  in  (5)  with  ”  0  to  get,  after  some  algebra. 


^W1'  ■  -  As#2’ '■ft  ■  <711 


where  we  have  only  retained  contributions  around  (0)  ,  k  )  .  Furthermore,  we 

o  o 

have  assumed  <0  /k  -  c  and  assumed  w  to  be  a  slowly  varying  function  of  z 
o  o  o  e 

in  the  sense  that 


2  2 
Id  \tf  /dz 
e 


« 


k  |dw  /dz| 
o  e 


(72) 


propogotlon 

direcfion 


Fig.  8.  Self-defocusing  of  a  beam  in  a  medium  with  a  negative  cubic 
nonlinearity.  The  beam  spreads  more  than  the  linear 
diffraction-limited  case. 
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_2  2  2  2  2 
In  (71),  denotes  the  transverse  Laplacian  8  /8x  +  8  /8y  .  (71)  with  P^  - 

0  is  identical  to  a  PDE  for  which  can  be  readily  solved  using  Fourier 

transform  techniques  £18]  to  yield  the  Fresnel  diffraction  formula  £4],  [19]. 

Thus,  in  (71) ,  the  first  term  on  the  RHS  is  due  to  diffraction,  while  the 

second  term  represents  the  nonlinear  contribution.  Eq.  (71)  has  the  same 

form  as  the  nonlinear  Schrodinger  equation  [2],  [4],  [17],  [20],  which  is 

used  to  explain  soliton  propagation  through  fibers  [21] . 

In  our  quest  for  the  expression  of  IY  |  that  doe3  not  depend  on  z,  we 

e 

may  substitute 

Ye  (x,  y,  z)  -  a  (x,  y)  exp  (-jkz)  (73) 

in  (71)  to  get 

vfa  -  -2ick  a  +  (P,k2/2)a3.  (74) 

t  o  oo 

Consider,  first,  the  case  where  we  have  one  transverse  direction,  viz., 
x.  Then  it  may  be  readily  verified  that  a  particular  solution  to  (74)  may 
be  expressed  in  the  form  [22] 

a(x)  «  A  sech  Kx  (75) 


where 


1/2 

A  -  <8K/P,k  )  ,  (76a) 

J  o 

K  -  l/(-2Kk  )1/2.  (76b) 

o 

We  note  that,  from  above,  K  <  0  and  p^  <  0  for  a  physical  solution.  Now  p^ 


35 


<  0  implies  >  0  (see  (12)),  which  is  in  agreement  with  our  heuristic 
description  for  self-focusing.  A  plot  of  the  envelope  is  shown  in  fig.  9. 

For  two  transverse  directions,  viz.,  x  and  y,  we  will  only  consider  the 
case  where  we  have  radial  symmetry,  and  express  the  transverse  Laplacian  in 
polar  coordinates: 

32/3x2  +  32/3y2  =  32/3p2  +  (1/p) 3/3p .  (77) 

Using  the  definitions 

a  =  (2K/P310 1/2  a 

p  -  (-l/2Kk  )  1/2fJ'  (78) 

eq.  (74)  may  be  recast  in  the  form 

32a/3p2  +  (i/p)3a/3p  -  a  +  a3  -  o.  (79) 

This  equation  ha3  no  analytic  solutions;  the  solutions  are  thus  obtained  by 
numerical  methods  [2],  (15],  [16]  and  are  shown  in  fig.  10.  These  are 

"multimodal"  in  nature,  the  mode  number  depending  on  the  initial  condition 
3(0)  . 

An  approximate  solution  for  an  arbitrary  initial  beam  profile  can  be 
found  by  starting  from  (71),  writing  down  the  "eikonal  equations"  during 
propagation  by  resolving  into  its  amplitude  and  phase,  and  assuming 
power  series  solutions  for  these  quantities  [23]  .  The  calculations  confirm 
periodic  behavior  if  we  consider  one  transverse  dimension;  but  are  outside 
the  scope  of  our  discussion.  We  end  this  section  by  showing  pictures  of  thi3 
periodic  behavior  for  an  initial  Gaussian  profile  (fig.  11a)  that  was 


Fig.  9.  Propagation  of  a  stationary  beam  (a)  having  the  form  as  given 
by  eqs.(75),  (76).  Fig.  (b)  was  generated  by  means  of  a 
split-step  algorithm  in  which  diffraction  was  accounted  for 
in  the  spatial  frequency  domain  and  nonlinearity  in  the 
space  domain.  The  stationary  beam  is  a  result  of  a  delicate 
balance  between  nonlinearity  and  diffraction. 


m  =  Mode  Number 


Fig.  10.  Numerical  solutions  to  eq.  (79)  showing  amplitude  profile 

for  higher-order  modes.  The  initial  values  for  the  respective 
modes  as  written  on  the  figure  ensure  that  the  solution  decays 
to  zero  at  infinity. 


Fig.  11.  (a)  Periodic  self-focusing  of  an  initial  Gaussian  profile  during 
propagation  through  a  self-focusing  medium;  (b)  mutual  attraction 
of  two  such  beams  in  aself-focusing  medium.  The  lower  picture  in 
(b)  is  an  equiphase  contour  to  highlight  the  mutual  attraction. 
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numerically  performed  by  employing  a  split-atep-type  angular  plane  wave 
spectrum  method  where  diffraction  is  accounted  for  in  the  spatial  frequency 
domain  (recall  our  comment  of  eq.  (71)  for  ”  0  before)  while  the 
nonlinear  effect  is  incorporated  in  the  space  domain  [24] .  The  interaction  of 
two  beams  adjacent  to  each  other  is  3hown  in  fig.  lib,  which  shows  mutual 
"attraction"  under  the  effect  of  nonlinearity. 


Optical  Bistabilit; 


Optical  bistability  refers  to  the  existence  of  two  stable  states  of  an 


optical  system  for  a  given  set  of  input  conditions.  It  is  interesting 
physically  because  it  represents  a  new  kind  of  nonlinear  system  in  optics. 
Such  an  effect  is  also  of  obvious  practical  interest,  as  it  offers  a  means 
of  realizing  all-optical  switching  for  optical  computers  [25] . 

The  typical  characteristics  of  an  optically  bistable  system  are 
illustrated  in  fig.  12.  A  gradual  increase  in  input  power  produces  a  3teady 
increase  in  the  output  power  or  intensity  until  reaching  a  critical  value 
where  the  output  jumps  up.  On  decreasing  the  input,  the  output  does  not 
immediately  fall  sharply  but  remains  on  the  upper  branch  of  the  curve  until 
the  input  is  reduced  to  a  lower  critical  value,  at  which  the  output  jumps 
down  again.  In  the  region  between  these  two  critical  points,  there  are  two 
stable  states  for  a  given  incident  power.  In  addition  to  being  hysteretic, 
the  system  also  shows  switching  transitions,  usually  at  both  edges  of  the 
bistable  region.  Also,  the  3tate  the  bistable  device  actually  assumes 
depends  on  the  direction  in  which  it  traverses  the  hysteretic  curve,  i.e.  it 
has  a  memory. 

In  general,  the  requirements  of  these  bistable  charcteristics  requires 
a  nonlinear  input/output  relationship  and  a  positive  feedback.  In  what 
follows,  we  will  first  set  forward  an  example  of  a  hybrid  bistable  system 
e.g.,  an  acoustooptic  bistable  device  to  point  out  the  requirements  for 
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0 


INPUT  POWER  - ► 


Fig.  12.  Typical  characteristics  of  an  optically  bistable  system. 

As  shown  in  the  sequence  of  pictures,  the  width  of  the 
bistable  region  may  be  changed  by  varying  the  feedback 
parameter. 
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bistability  mentioned  above.  This  will  be  followed  by  examples  of  a 
nonlinear  Fabri-Perot  cavity  and  a  linear  nondispersive/nonlinear 
(dispersive)  interface-type  arrangement,  both  of  which  exhibit  optical 
bistability.  An  excellent  reference  on  optical  bistability  is  the  book  by 
Gibbs  [26]. 

5 . 1  Acoustooptic  Bistability 

If  the  diffracted  light  from  an  acoustooptic  device  operating  in  the 

Bragg  regime  is  detected,  amplified  and  fed  back  to  the  transducer  driving 

the  acoustooptic  cell,  a  bistable  device  results.  Fig  13  shows  the 

experimental  arrangement  [27]  .  An  acoustooptic  cell  is  driven  with  a  40  Mhz 

generator  which  creates  the  acoustic  grating.  The  first-order  diffracted 

light  from  a  He-Ne  laser  is  detected  and  fed  back  to  the  transducer  in  phase 

with  the  external  acoustic  "input".  If  we  denote  the  incident  light  field 

by  y.  ,  the  undiffracted  and  diffracted  light  fields  by  \ir  and  w 
xnc  o  1 

respectively,  then  the  "output"  Y1  is  given  by 

*1  -  Vinc  sin  (a/2),  (80) 

where  a  is  the  peak  phase  delay  encountered  by  the  light  in  the  interaction 
region  and  is  proportional  to  the  amplitude  S  of  the  sound  field.  Note  that 
Vn  “  y.  cos  (a/2),  to  maintain  conservation  of  energy.  Note  that  the 
assumption  made  in  deriving  (8))  is  that  the  sound  pressure  («a)  remains 
constant  during  the  interaction. 

If,  now,  the  first  order  diffracted  light  intensity  is  detected  by  a 
photodetector  and  the  resulting  electrical  signal  amplified  (by  passing  it 
through  an  amplifier  with  a  gain  constant  |$)  and  fed  back  to  the  transducer 
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Fig.  13.  (a)  Acoustooptic  bistable  device  operating  in  the  Bragg  regime; 
(b)  hybrid  bistable  device 
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in  unison  with  external  bias  (a  ) ,  we  have  a  system  with  a  nonlinear  input 

o, 

(a)  -  output  (y^)  relationship,  and  having  a  positive  feedback.  The 
effective  a  scattering  the  light  in  the  acoustooptic  cell  may  be  written  as 

a  -  &o  +  (5  ly^J2,  (81) 

where  j5  denotes  the  product  of  the  gain  constant  (3  of  the  amplifier  and  the 

conversion  efficiency  of  the  photodetector. 

Note  that  under  the  feedback  action,  a  can  no  longer,  in  generality,  be 

treated  as  a  constant.  In  fact  a  can  be  treated  as  constant  during 

interaction  iff  the  interaction  time,  given  as  the  ratio  of  the  laser  beam 

width  (r  1  mm)  and  the  bulk  speed  of  sound  in  the  cell  (-  2000  m/3),  is  much 

smaller  compared  to  the  delays  incorportated  by  the  finite  response  time  of 

the  photodetector,  the  r.f.  sound  cell  driver  and  the  feedback  amplifier,  or 

any  other  delay  line  that  may  be  purposely  installed  (e.g.,  an  optical  fiber 

or  coaxial  cable)  in  the  feedback  path.  We  will  consider  thi3  case  only. 

Based  on  the  discussion  above,  we  conclude  that  corresponding  to  a 

2  2 

fixed  input  (<*  a0) ,  the  output  ^(or  |y  I  *  y^  xf  we  assume  y^  to  be  real) 

will  undergo  a  series  of  iterations  at  every  instant  a  is  updated  through 

2 

the  feedback  action.  The  value  of  y^  after  n  iterations  may  be  written, 
using  (80)  and  (81),  as  [28] 

V,  (n+l)  -  y.  sin  ((a  +  jSy.  (n) ) /2]  .  (82) 

1  Txnc  o  1 


The  analysis  of  (82)  involves  a  knowledge  of  the  theory  of  autonomous 
dynamical  systems  of  which  only  rather  recent  analysis  has  been  performed 
for  quadratic  maps  of  the  type 
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xn+l  "  ^nU‘Xn>' 


(83) 


and  is  out  of  the  scope  of  our  discussion.  Suffice  to  state  that  the 
behavior  of  (83)  is  dynamically  simple  for  0  <  (l  <  3  but  chaotic  when  p.  >  4 
[28],  [29].  In  fact,  the  map  undergoes  a  series  of  period-doubling 
bifurcations  en  route  to  chaos  as  the  value  of  p.  increases. 

In  what  follows,  we  will  discuss  the  steady-state  behavior  of  (82)  (n  -* 

«>)  for  different  values  of  .  Rather  than  go  through  an  extensive  stability 
analysis  of  the  steady  state,  we  will  draw  the  analogy  between  the  operation 
of  the  bistable  device  with  an  electronic  Schmitt  trigger  [30]  which  also  exhibits 
hysteretic  properties,  and  give  an  equivalent  circuit  representation.  First, 
note  from  (82)  that  the  steady  state  relationship  reads 


V1  “  VincSin2  [(“°  +  PVi)/2] 


(84) 


and  is  plotted  in  fig.  14  for  different  values  of  {5.  Note  that  (5  *  2.2  is 
the  lowest  admissible  value  of  the  feedback  parameter  to  observe  any 
hysteresis,  with  increase  in  (5,  the  area  under  the  hysteresis  curve 
increases,  and,  more  noticeably,  the  lower  threshold  moves  to  the  left.  This 
makes  sense  since,  upon  differentiating  (84)  w.r.t.  a^,  we  get 


dy.  /da  -  (y.  /2)  (sin  a) / [1- (p/2)  (sin  a)]; 

x  u  inc 


(85) 


hence  dw,  /da  will  become  infinite  when 
1  o 


<{5/2)  sin  o-l. 


(86) 


Fig.  14.  Sketch  of  the  steady-state  relationship  for  the  hybrid 
acoustooptic  bistable  device.  Transitions  are  marked  by 
dotted  lines  and  path  abed  completes  the  hysteresis  cycle 
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which  implies  ^  >  2  for  the  onset  of  hysteresis. 

In  passing,  we  remark  that  in  a  practical  Bragg  operation,  seme  amount 

of  light  is  also  diffracted  into  higher  orders,  noticeably  the  second  order. 

Issues  such  as  improvement  of  performance  by  reducing  the  second-order 

diffraction  or  the  effect  of  a  second-order  feedback  are  current  topics  of 

interest.  We  also  mention  that  similar  hy3teretic  behavior  may  be  observed 
2  - 

by  treating  V^nc  and  P  as  the  inputs  and  treating  the  other  variables  a3 
parameters  [27] . 

For  our  configuration,  the  equivalent  circuit  in  terms  of  an  electronic 
Schmitt  trigger  is  shown  in  fig.  15  [28].  Readers  familiar  with  an  OP-AMP 
realization  of  the  Schmitt  trigger  may  recognize  that  the  transfer 
characteristic  of  the  circuit  when  the  loop  gain  is  greater  than  1  has  a 
3hape  (see  fig.  16)  similar  to  the  curves  in  fig.  14.  The  part  of  the 
curve  with  a  negative  slope  is  unstable,  and  there  are  two  stable  3tate3. 
Transition  from  one  state  to  the  other  occurs  at  the  upper  and  lower 
thresholds  of  the  input.  The  optically  bistable  device  operates  in  a 
similar  way;  with  transitions  along  the  dotted  lines  in  fig.  14. 

5.2  The  Nonlinear  Fabri-Perot  Cavity 

Consider  a  CW  laser  beam  injected  into  an  optical  cavity  which  is  tuned 
or  nearly  tuned  to  the  incident  light.  In  general,  the  incident  field  is 
partly  transmitted,  partly  reflected  and  partly  absorbed.  When  the  cavity 
contains  an  absorbing  material  resonant  or  nearly  resonant  with  the  incident 
electric  field,  the  transmitted  power  becomes  a  nonlinear  function  of  the 
incident  power.  The  behavior  of  the  system  is  determined  by  a  parameter 
which  depends  on  the  unsaturated  absorption  coefficient  of  the  sample  per 
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3  =  ^i nc/2)sin  ^ 


Fig.  15.  Comparison  of  the  hybrid  acoustooptic  device  with  an  electronic 
Schmitt  trigger.  The  small  signal  gain  with  feedback  is  given 
by  eq.  (85).  Unlike  the  electronic  Schmitt  trigger,  the  gain 
characteristics  of  the  amplifier  do  not  possess  a  hard 
nonlinearity. 


Fig.  16.  Transfer  characteristics  of  an  electronic  Schmitt  trigger  using 
the  model  of  an  ideal  OP-AMP.  V  is  the  output  while  V. 
represents  the  input,  and  should°be  compared  to  -a  in  £ig.  15. 
Comparison  with  fig.  14  reveals  the  essential  differences: 

(a)  the  curve  seems  reflected  about  a  vertical  line,  this  is  due 
to  the  fact  that  the  inputs  are  similar  upto  a  negative  sign, 
and  (b)  the  upper  and  lower  steady  states  are  well-defined,  due  to 
the  hard-limiter  model  of  an  ideal  OP-AMP. 
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unit  length,  its  length  and  on  the  mirror  transmittance.  If  one  increases 
the  value  of  this  parameter  above  a  certain  threshold,  the  steady-state 
input-output  curve  becomes  S-shaped,  indicating  bistability  and  hysteresis. 
The  threshold  value  depends  on  the  cavity  mistuning,  the  atomic  detuning, 
the  inhomogeneous  linewidth  and  the  type  of  cavity,  to  name  a  few.  When  the 
incident  field  is  in  perfect  resonance  with  the  atomic  line  so  that 
dispersion  does  not  play  a  role,  we  have  purely  absorptive  bistability.  On 
the  other  hand,  if  the  atomic  detuning  is  so  large  so  that  absorption 
becomes  negligible,  we  have  purely  dispersive  bistability  [31] . 

In  order  to  describe  theoretically  the  phenomenon  of  bistability,  we 
will  consider  a  unidirectional  ring  cavity  as  shown  in  fig.  17.  For 
simplicity,  we  assume  that  mirrors  3  and  4  have  100%  reflectivity.  We  will 
denote  the  reflection  and  transmission  coefficients  of  mirrors  1  and  2  a3  R 
and  T  respectively. 

Now,  the  dynamics  of  the  "active"  cavity  between  mirrors  1  and  2  may  be 
modelled  in  terms  of  the  Maxwell-Bloch  equations  [26] : 

3E  /dt  +  c  8e  /8z  -  -gP  ,  (87a) 

e  o  e  e 

3p  /3t  -  (pt/R )  E  5-  jAP  ,  (87b) 

e  e  e 

3D/8t  -  (JI/2R)  (E  P*  +  E*P)  -  y  (D-N/2 )  ,  (87c) 

e  e  e  1 

where  E  and  P  are  the  slowly  varying  envelopes  of  the  electric  field  and 
e  e 

macroscopic  polarization,  cq  represents  the  linear  phase  velocity,  R  -  h/2K 
where  h  is  Planck's  constant  and  )I  is  the  modulus  of  the  dipole  moment  of  the 
atoms.  Furthermore,  D  is  one  half  the  difference  between  the  populations  of 


Fig.  17.  The  nonlinear  Fabri-Perot  modelled  by  a  ring  cavity.  The 
nonlinear  meterial  is  between  mirrors  1  and  2.  E.,  E  and 

E  denote  Che  incident,  reflected  and  transmitted  electric 
fields  as  shown. 
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the  lower  and  upper  level,  g  is  a  coupling  constant  proportional  to  the 

frequency  (0  of  the  incident  field,  A  -  (a )  -(0  )-jy_  where  o  is  the 
o  a  o  2  a 

transition  frequency  of  the  atoms  and  y  ,  y2  are  the  inverses  of  the  atomic 
relaxation  times  and  respectively. 

In  what  follows,  we  will  analyze  the  absorptive  and  dispersive  cases 
separately.  For  the  first,  we  will  bring  out  the  steady  state  relationship 
and  draw  the  hysteresis  curve  to  bring  out  the  similarity  (and  differences) 
with  the  case  considered  in  Section  5.1.  For  the  second,  we  will  get  into 
a  more  detailed  analysis  to  bring  out  the  effects  of  feedback  gain  and  the 
onset  of  instabilities. 

A.  Absorptive  Case 

Here  (0  »  (0  so  that  jA  (see  (87b))  =  y,,;  E  ,  P  may  be  considered  real 

a  o  zee 

without  loss  of  generality  to  simplify  (87c) .  The  coherent  CW  field  E  ^ 

enters  into  the  cavity  from  the  left.  The  cavity  imposes  the  two  following 

relations  between  E  ,,  the  transmitted  field  E  and  the  fields  E  (0,t)  and 

ei  et  e 

E  (L, t)  where  L  is  the  length  of  the  atomic  sample: 
e 

E  .  (t)  =  T1/2E  (L,t) , 
et  e 

1/2 

E  <0,t)  -  T  E.  +  RE  (L,t-At).  (88) 

e  i  e 

In  (88),  At  •  (2£+L)/c  is  the  time  light  takes  to  travel  from  mirror  2  to 

o 

mirror  1  and 


R  +  T 


1. 


(89) 
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The  second  equation  in  (88)  is  a  boundary  condition  characteristic  of  the 
ring  cavity  where  the  second  term  on  the  RHS  describes  a  feedback  mechanism 
so  essential  for  bistability. 

Consider  the  steady  state  3e  /3t  -  3P  /3t  -  3D/3t  -  0.  The  set  of 

0  0 

equations  (87)  may  be  readily  simplified  to  yield 


dE  /dz 
e 


2 

-(HgN/2Rc  Y  )  E  /  [  1+E2/  (H  y.y^/\L)]. 
o  2  e  e  12 


(90) 


Eq.  (90)  shows  how  the  "effective"  propagation  constant  for  the  electric 
field,  and  hence,  the  refractive  index,  depends  on  the  electric  field 
itself,  and  thus  brings  out  the  physics  of  the  nonlinearity  constants  viz. 
P  ,  ,  n^  and  y  introduced  in  Section  2. 

For  notational  convenience,  we  introduce  the  parameter 


«c  *  ngN/2Rc0Y2 

and  a  normalized  field 


(91) 


>lEe/R(YlY2) 


1/2 


(92) 


to  recast  (90)  into 

dy  /dz  -  -  a  y  /(l+y2),  (93) 

0  C  ©  0 

which  can  be  easily  integrated  and  evaluated  at  z-L.  The  result  is 

4n  (y  (0)  /\jf  (L) )  +  (1/2)  (0)  -  vfi2(L))  -  a  L. 

e  e  e  e  c 


(94) 
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Now,  we  define  normalized  incident  and  transmitted  amplitudes  as 


V, 


ei,  t 


HE  ,  JK11, Y?T) 
ex,  t  j.  ^ 


1/2 


(95) 

I 


Then  (88)  in  the  steady  state  transforms  to 


et 


Ve<L) , 


V0) 


Ty  .  +  Ry 


ei 


et 


(96) 


With  (96),  eq.  (94)  may  be  rewritten  as 


«n[l+T(Vei/Vet 


1)]  +  (vj/2){(l  +  -  DJ2  -1) 


a  L. 
c 


(97) 


In  the  so-called  double  mean-field  limit  [31] 


a  L  «  1,  T  «  1  with  CXL/2T  -  C  -  constant,  (98) 

c 

(97)  has  the  form 

Yei  "  \t  +  2C*et/(1+  V2)  (99) 

and  is  plotted  in  fig.  18.  Note  that  the  graph  is  an  S-shaped  curve 
(similar  to  the  acoustooptic  case)  which  leads  to  a  hysteretic  cycle. 

We  remark  that  the  bistability  arises  from  the  combined  action  of  the 
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nonlinear  input-output  relationship  of  the  sample  medium  (see  (94))  and  the 
feedback  of  the  mirrors.  This  feedback  action  is  essential  since  it  is  easy 
to  show  that  a  plot  of  (97)  for  T  *  1  (no  reflection)  does  not  possess  the 
characteristic  S-shape  any  more. 

Now,  introducing  perturbations  of  the  steady-state  solution  to  the 
Maxwell-Bloch  equations,  one  can  perform  a  rigorous  (although  local) 
stability  analysis  for  the  resulting  linearized  system.  Stability  then 
requires  that  the  perturbations  decay  exponentially  in  time;  hence  the 
eigenvalues  of  the  linearized  problem  should  have  nonpositive  real  part3.  A 
detailed  analysis,  which  is  done  in  [31],  is  out  of  the  scope  of  this 
discussion. 

B.  Dispersive  (Kerr)  Case 

For  an  analysis  of  Kerr-type  bistability,  we  return  to  the  Maxwell- 
Bloch  equations  (87)  and  assume  CO  -  0)  »  Y-  For  a  simplistic  approach 

cl  O  ^  • 

and  to  decouple  the  equations  in  the  easiest  possible  way,  we  can  eliminate 

D  between  (87b)  and  (87c)  at  around  the  steady  3tates  for  P  and  D  and 

e 

substitute  in  (87a)  to  obtain 

3y  /3t  +  c  /3z  -  j(gnN/R|A|)y  -  j  (gja.3^  N/2y  H3  |  A I  3 )  \|/  2y*  (100) 

0  0  6  ©  fc  1  6  c 

where  we  have  written  y  for  E  . 

e  e 

It  can  be  verified,  after  some  algebra,  that  if  yg  satisfies  (100),  then 
the  total  wavefunction  y  satisfies 

32y/3t2  -  co2a2y/3z2  -  A^  +  (2^73)3  V/3t2,  (101) 
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where 


A.  -  -gpNk  c  /K|A|, 
l  o  o 


(102a) 


P3  -  -4g)i3Y2N/71fi3|A|3  kQCQ, 


(102b) 


if  we  restrict  ourselves  around  an  operating  frequency  CO  .  Note  (102b)  shows 

o 

the  physics  of  the  origin  of  the  Kerr-type  nonlinearity.  Also,  an  inherent 
dispersion  appears,  as  depicted  through  A^ .  In  fact,  the  linear  dispersion 
relation  can  be  obtained  by  suppressing  the  nonlinear  term  and  substituting  w 
-  exp  j  (cot  -kz)  ,  and  reads  [2] 


2 

Eq.  (103)  resembles  a  waveguide-like  dispersion  relation  if  A  =  -co  ,  where 

1  c 

co  is  the  (linear)  cutoff  frequency  of  the  waveguide.  We  have  thus  elegantly 
c 

reduced  the  complicated  Maxwell-Bloch  equations  to  a  form  similar  to  (4b) 
(which  was  heuristically  derived  from  a  simple  physical  picture  of  phase 
velocity  dependence  on  amplitudes) . 

We  will  u3e  (101)  as  a  model  to  gain  insight  into  the  instabilities 
encountered  in  the  Fabri-Perot  cavity  with  feedback,  as  well  as  to  explain 
hysteresis  and  bistability  for  the  interface  problem  in  the  next  section. 

In  what  follows,  we  provide  a  simple  analysis  of  instantaneous 
optical  bistability  based  on  the  model  of  the  nonlinear  medium  as  described 
above  by  treating  it  as  a  degenerate  parametric  amplifier  of  any  noise 
present  in  the  cavity.  Depending  upon  the  phase  of  the  noise  component,  the 
net  parametric  round  trip  gain  may  be  real  and  exceed  unity  magnitude. 
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Under  such  conditions,  we  expect  instabilities. 

From  a  physical  point  of  view,  the  nonlinear  medium  acts  as  a  degenerate 
traveling  wave  parametric  amplifier  for  any  noise  in  the  system  where  the 
"pump",  which  acts  as  frequency  2C0q/ exists  virtually  through  the  action  of 
the  cubic  nonlinearity  on  the  signal  itself  [32]  .  In  our  quest  for 
instabilities,  we  investigate  if  it  is  possible  for  a  noise  component  of  a 
certain  phase  to  exist  such  that  after  parametric  amplification  in  the  medium 
and  completion  of  its  round  trip  through  the  ring  cavity,  this  component  will 
return  in  the  nonlinear  medium  in  the  same  phase  as  when  it  started.  In 
cases  where  this  possible,  we  would  then  expect  instability  if  the  roundtrip 
gain  exceeds  unity. 

An  alternate  way  of  writing  (100)  or  (101)  is  to  express  y  as  a  phasor: 

w  =*  Re  {y  exp  j(0  t) ;  =  \U  (z),  (104) 

p  o  J  p  p 

where  we  have  assumed  CW  propagation.  Substituting  this  in  (101)  yields 

d2yp/dz2+  <k2  +  A1/cQ)Vp  -  (P3ko/2)VpVp  (105> 

if  we  restrict  ourselves  around  (0  and  denote  k  ■  to  /c  .  Now,  the  total 

o  o  o  o 

wavefunction  yp  may  be  visualized  as  the  sum  of  a  signal  and  noise  at  the 
same  frequency  and  decomposed  as 

Vp  -  Vpl  +  «Vp2,  U06) 

where  e  is  a  small  quantity.  An  order-by-order  analysis  of  (105)  after 
substitution  of  (106)  gives 
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e°:  d2ypl/dz2  +  l^+A1/c2)Vpl  -  <IVp/2>Vp2V 


E1:  d2yp2/dz2  +  l^+Aj/c2)^  -  (P3k2/2)  ivp2yp2%2yplyp*vp2) 


Assuming  an  undepleted  pump  \i/  ,  we  next  write,  in  (107a), 

pi 


Vpl  “  Vel(2)  ®Xp  (_jkoz) 


to  determine  the  modified  propagation  constant  k^ 


k'2  -  (k2  +  A./c2)  -  (P,k2/2) IV  |2. 
o  o  lo  *3o  Tel 


Note  that  the  effective  propagation  constant  is  modified  from  kQ 
contributions  from  the  dispersion  (A^)  and  the  nonlinearity  (P^) 
study  the  behavior  of  V  o  (the  "noise"  component),  we  3et 

p2 

VP2  ‘  Ve2U)  exp 

and  substitute  in  (107b)  to  obtain 


d»e2/dz  -  j(P3k2/Ek')(veiye2tye*yp2lyel, 

where  we  have  used  (109)  and  the  "slowly  varying"  assumption 

|d2y  /dz2|/|dv  ,/dzl  «  k' 

e2  ez  o 


(107a) 

(107b) 

(108) 


(109) 

due  to 
Now  to 

(110) 


(111) 


(112) 
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to  simplify. 

In  order  to  track  the  amplitude  and  phase  of  the  "noise"  component 

separately,  we  write 


V, 


el 


V 


Ve2  -  a2 (z) exp  -j$2  (z) , 


(113) 


where  we  have  assumed  the  pump  amplitude  to  be  real  without  loss  of 
generality.  Substituting  (113)  in  (111),  we  get  the  pair  of  coupled 
equations 


da2 / dz  -  -(P3k^a^/4kMa2  sin  2$2, 
d<|>2/dz  -  -(P3k^a^/4kM  (l+cos  2<J»2>  , 


(114a) 


(114b) 


upon  separating  real  and  imaginary  parts.  Now  (114b)  may  be  readily 
integrated  to  give 


<t>2  (z) 


tan  1  [(P  kV/2k')z], 

j  o  i  o 


(115) 


where  we  have  assumed  $2(0)  “  0  for  simplicity,  and  where  k^  is  defined 
through  (109)  and  (113) .  Substitution  of  (115)  in  (113)  and  integration 
yields,  after  some  algebra, 


22  22  1/2 

a,(z)  -  a  (0)  [l+(P,kV/2k')  *  ]  . 

2  2  J  O  I  O 


(116) 
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With  (115)  and  (116),  we  can  write 


rp2 


-  iyp2(0) i (i+(az)2]1/2exp 


j (tan  -  k' z) , 
o 


(117) 


where 

« -  P3koi^i|/2it;-  'n8» 

Eq.  (117)  shows  that  the  "noise"  component  is  amplified  as  it  travels 
through  the  nonlinear  medium  of  length  L.  In  fact, 

2  1/2  -1 

V  ,(L)  -  |y  -<0) I [l+(OL)  ]  exp  -j (k'L-tan  OL) ,  (119) 

pZ  pZ  O 

so  that  the  "noise"  component  phasor  after  one  round  trip  read3 

2  1/2 

y  (2L+2t)  -  R|y  (0) | [l+(OL)  ]  exp  -j(k'L+k  (L+2f)-tan  OL) (120) 
pz  pz  o  o 

where  we  have  taken  into  account  the  phase  change  during  propagation  through 
the  path  2341  (see  fig.  17)  and  introduced  the  effects  of  reflection  by 
mirrors  2  and  1 . 

The  conditions  for  oscillation,  or  instability,  (like  in  the 
acoustooptic  case)  demand  that  (a)  the  total  round  trip  phase  change  must 
be  an  integral  multiple  of  2jc  (positive  feedback)  ,  and  (b)  the  round  trip 
(real)  gain  A  must  exceed  1.  The  first  condition,  in  our  context,  reads 

(k'L  -  tan  ^QL)  +  k  <L+2£)  -  2rwt,  m  integer. 


(121) 
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while  the  second  condition  yields 

2  1/2 

A  -  R[l+  (QIi)  ]  '  >  1,  (122a) 

implying 

1  +  (OL)2  >  ( 1/R) 2 .  (122b) 

Using  (121),  the  gain  condition  (122b)  may  be  rephrased  as 

R  >  cos[(k'-k  )  L  +  2L  (L+t)  -  2imt]  .  (123) 

o  o  o 

The  second  term  in  the  argument  of  the  cosine  in  (123)  denotes  the  total 
phase  shift  per  round  trip  in  the  Fabri-Perot  in  the  absence  of  the  nonlinear 
sample;  while  the  first  term  denotes  the  extra  phase  3hift  due  to  dispersion 
(A^)  and  nonlinearity  (P^)  in  the  sample.  (121)  and  (123)  should  serve  as 
good  "design"  equations  in  the  construction  of  a  Kerr-type  nonlinear  Fabri- 
Perot  exhibiting  bistability.  Our  simplified  analysis  is  in  good  agreement 
with  rigorous  theories  predicting  bistability  in  such  a  device  [33] . 

5 . 3  The  Linear/Nonlinear  Interface 

Consider  the  semi-infinite  linear  and  nonlinear  regions  as  shown  in  fig. 
19  with  a  common  interface  at  z  -  0.  Assume,  also,  a  plane  wave  incident 
from  the  linear  medium  to  the  nonlinear  medium  at  "grazing  incidence."  Under 
the  assumptions  that  the  reflected  wave  is  a  plane  wave  and  that  the 
transmitted  field  is  either  a  propagating  plane  wave  or  an  evanescent  field 
which  decreases  monotonically  with  distance  from  the  interface  z,  but  is 
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Fig.  19.  Wave  diagram  of  the  nonlinear  interface  and  for  "grazing" 
incidence.  The  ray  traces  in  the  nonlinear  medium  (below 
the  interface)  represent  (a)  the  traveling  plane  wave  (PW), 

(b)  the  longitudinally  inhomogeneous  traveling  wave  (LITW) 
and  (c)  the  totally  internally  reflected  wave  (TIR) .  The 
nonlinear  material  is  assumed  to  possess  Kerr-type  nonlinearity. 
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independent  of  x  and  y  coordinates,  the  scalar  wave  equation  can  be 
analytically  solved  in  closed  form  and  satisfies  all  the  relevant  boundary 
conditions.  Solutions  can  be  obtained  for  both  positive  and  negative  values 
of  Pj*  The  theory  predicts  hysteresis  and  bistability  for  the  reflection 
coefficient  when  plotted  as  a  function  of  the  input  amplitude,  and  is  drawn 
in  fig.  20  [34]  .  Later,  numerical  simulations  were  done  to  analyze  the  reflection 
of  a  Gaussian  beam  from  the  interface,  also  at  grazing  incidence  [35],  and  some  of 
the  predicted  effects  were  experimentally  observed  [36] . 

In  what  follows,  we  will  consider  a  simpler  model  where  a  plane  wave  is 
normally  incident  at  the  interface  between  a  linear  nondispersive  medium  and 
a  nonlinear  dispersive  medium  (see  fig.  21)  [37],  Normal  incidence  makes  the 

computations  simpler,  and  dispersion  is  required  to  ensure  bistable 
operation.  Dispersion  may  be  provided  by  having  a  waveguide-type  structure 
containing  the  nonlinear  medium.  We  could  also  think  of  the  nonlinear  medium 
to  be  excited  at  a  frequency  much  different  from  the  transition  frequency  of 
the  atoms  so  that  the  material  appears  dispersive  to  the  incoming  field  (see 
section  5.2). 

Before  proceeding  with  the  analysis,  let  u3  remark  that  the  linear 
region  will  be  modelled  by  a  linear  waveeqution  of  the  type 

d2y/3t2  -  cq2  92y/8z2  -  0  (124) 

and  the  nonlinear  dispersion  region  by  an  equation  similar  to  (101) : 

a2y/at2  -  d2  a2y/az2  -  A  y  +  <2p3/3)aV/at2.  (125) 

We  have  replaced  the  linear  phase  velocity  by  in  the  nonlinear  medium  to 
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Fig.  20.  Dependence  of  the  reflectivity  r  on  the  incident  light  intensity 
lE^I  at  different  glancing  angles  in  the  case  of  negative 
linear  mismatch  of  susceptibilities  and  n^  greater  than  zero. 
Increasing  label  numbers  refer  to  decreasing  ratios  of  the 
glancing  angle  to  the  critical  glancing  angle,  which  is 
determined  from  the  linear  mismatch  of  the  susceptibilities. 
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bring  out  the  effect (a)  of  linear  phase  velocity  mismatch.  A  time-harmonic 
incident  wave  will  be  assumed,  and  in  the  propagation  mode,  the  wave 
amplitude  in  the  nonlinear  medium  will  be  determined  through  the  amplitude 
transmission  coefficient 

T  -  2TJ  /  (Tj  +  TV)  =  2k '  /  (k  +k '  )  ,  (126) 

o  o  o  o  o  o 

where  kQ,  k^  refer  to  the  propagation  constants  in  the  linear  nondispersive 
and  nonlinear  dispersive  regions  respectively  [38]  .  Sometimes,  there  is  no 
propagation  in  the  nonlinear  region;  then  only  evanescent  modes  can  exist. 
The  wave  amplitude,  as  will  be  shown  below,  decays  with  distance;  however, 
this  decay  is  not  exponential.  First  hand  derivations  will  be  done  for  this 
case.  The  output  is  a33umed  to  be  the  amplitude  of  the  electric  field 
beyond  a  certain  specific  distance  into  the  nonlinear  medium.  The  steady 
3tate  input  (incident  field) -output  characteristics  will  be  plotted  and 
their  stability  studied  to  bring  out  the  hysteretic  and  bistable  effects. 

Starting,  then,  from  (125),  and  assuming  a  phasor  form  for  y  as  before 
(see  (103)),  we  obtain,  analogous  to  (104), 

d2y  /dz2  +  k2y  +  (A  /2q)  \|/  -  (3,£2/2)\(/2\|/*  =  0,  k  -  0)  /c  .  (127) 

p  o  p  1  °  p  r3  o  P  P  o  o  o 

A  particular  traveling  wave  solution  of  the  above  has  the  form 

\y^(z)  -  Afc  exp  -jk^z,  Afc  constant  (128) 


with 
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k  =  ±  (1/2  )  f(02  -  (o2 
t  c  o  c 


<P  C02/2)  IA  I2]1'  2 
3  O  t 


(129) 


In  (128),  (129)  the  subscript  ’t1  denotes  transmission.  Also,  in  (129),  A^ 

2 

ha3  been  replaced  by  -co^  in  accordance  with  the  discussion  following 

sq.  (103).  The  transmission  coefficient  may  now  be  found,  using  (126),  (129), 

and  the  fact  A^  -  A,t  where  A.  is  the  incident  wave  amplitude  (assumed  real) 

tii 

from  the  linear  nondispersive  region,  in  a  trascendental  form  as 


X  =  2/ [1  ±(k  /k  )  {1 
o  o  L 


(CO  /CO  ) 
c  o 


(P3/2)A2t2}172]  . 


(130) 


The  presence  of  nonlinearity  changes  the  cutoff  frequency  of  the 
waveguide.  To  see  this,  note  that  for  the  condition  of  propagation  to  be 
satisfied  the  quantity  in  braced  brackets  in  (130)  must  be  greater  than 
zero;  and  may  be  achieved  with  p^  <  0  even  if  the  linear  cutoff  frequency 
(CDJ  is  greater  than  the  operating  frequency  «D  )  .  Under  these  conditions, 
we  may  3ay  that  a  strong  enough  incident  wave  can  "push"  through  a  waveguide 
in  linear  cutoff.  Alternatively,  and  in  the  case  of  our  interest,  there  is  no 
propagation  (kt  imaginary)  for  P^  >  0  (n^  <  0)  even  for  coc  <  0)q  if  the 
incident  field  is  strong  enough  such  that 


P  A2/2  >  [1  -  (co  /CO  )  2 ]  /X 2  . 
3  i  co 


(131) 


In  this  ca3e  it  is  easy  to  see  that  the  "effective"  cutoff  frequency  is 
increased  to  a  value  greater  than  the  operating  frequency. 

While  it  is  easier  to  predict  the  amplitude  of  a  propagating  wave  in 
the  nonlinear  medium  (using  (130) ;  ;he  evanescent  solutions  in  the 
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nonpropagating  case  has  to  be  determined  from  first  principles.  Since  the 
propagating  constant  is  now  imaginary,  we  may  write,  from  (129), 


-  2  2  1/2 
kt(z)  -  -jkQ(33  !At  (z)  I  -  {l-(Oc/03o)  }] 


»  - jS  (z)  . 


(132) 


Since  A  decays  in  the  nonlinear  dispersive  medium,  we  may  write  A  (z+Az)  = 
t  ^ 

A^tzjexp  -5(z)Az,  which,  in  the  limit  Az  ->  0,  read3 


dAt(z)/dz  =  -5 ( z ) A^ ( z ) 


At  (z) 


-y]1/2At(z) 


(133a) 


where 


P  £  P,/2;  Y  =  l-(0)  /to  )2.  (133b) 

r  r3  co 

[It  may  be  formally  shown  that  P,  y  >  0  to  achieve  bistability  [37);  however, 
this  i3  out  of  the  scope  of  this  rather  elementary  discussion] . 

To  solve  for  At<z),  we  first  resolve  it  into  its  amplitude  and  phase  by 
writing 


At  (z)  -  At  (z)  exp  j<j>  (z)  . 


(134) 


Substitution  in  (133)  and  some  algebra  [37]  shows  that 
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$  ( z )  -0(0)  -  constant  ; 


At(z)  -  (Y/P)1/2/sin  71/2{koz  +  (l/Y1/2)sin"1((Y/p)1/2/At(0))] 


(135) 


A  plot  of  the  solution  A^(z)  is  shown  in  fig.  22.  Note  that  A^(z)  takes  on 
its  maximum  value,  viz.  A^CO),  at  z=0,  and  decays  (nonexponentially)  to  a 


critical  "threshold"  value  (Y/(J)  at 


-1  1/2  —  1/2- 
z  =  z  =  [rt/ 2  -  sin  (  (7 /(3 )  /A  (0))]/Y  k  . 
o  to 


(136) 


1/2  - 

Incidentally,  ( Y /P )  Is  the  minimum  value  of  A^(z)  for  nonpropagation.  It 

can  be  readily  checked  that 


y(z,t)  »  (Y/P>  exp  j<t>Q  exp  j(0Qt 


(137) 


is  a  solution  of  (127)  for  z  >  z  and  the  k  =  0  for  this  case.  Thus,  A  (z) 

o  t  t 

decays  first  as  in  (135)  till  it  reaches  the  critical  value  from  which  point 
it  remains  constant  (see  fig.  22) . 

A  discussion  of  boundary  conditions  at  the  interface  (z=0)  for  the 
nonpropagating  case  is  now  in  order.  This  can  be  done  rigorously  by 
imposing  the  conditions  that  the  wavef unctions  and  their  spatial  derivatives 
are  continuous  across  the  boundary.  Alternatively,  this  can  be  done  by 
realizing  that  eqns.  (126),  (130)  hold  at  z»0+.  This  gives 


—  —  ?  2  —  ?  i/2  1/2 

A  (0)  -  t  [-{  (k  /k  )  -y}  +  {<<k  /k  )  -Y>  +  1 6  A^  P  ( k  /  k  )  }  1  /  j  /  2  J1  ]  , 

t  OO  OO  iOO 


(138) 


Fig.  22.  Plot  of  the  field  amplitude  into  the  nonlinear  dispersive  medium 
in  the  nonpropagation  mode  (see  eq.  (135)). 


71 


tan  <j>(0)  -  C  C-{  /k  )2+y}  +  {  <  (k  /k  )2-y>2  +  16a2  P(k  /k  ) 2} 1U} /20] 1/2 . 
o  o  oo  100 

(139) 

The  reflection  coefficient  has  a  magnitude  1,  which  is  not  surprising. 

Now,  assume  the  output  plane  to  be  defined  as  z=z^  in  the  nonlinear 
medium.  Fig.  23  shows  the  plot  of  A^_(z^)  vs .  A,  for  a  set  of  suitably 
chosen  parameters  P,  y,  k  ,  k^  and  z^.  It  comprises  two  main  parts:  the 
lower  part  (I  and  II)  which  represents  the  steady-state  solution  for  the 
transmission  case,  while  the  upper  part  (III  and  IV)  represents  the  solution 
for  the  nonpropagation  case.  Note  that  in  the  nonpropagation  mode,  A^(z^) 
remains  constant  with  increase  in  A.  the  corresponding  z  (which  depends  on 
A^  ( 0 ) ,  and  hence  on  AJ  is  less  than  z^.  The  figure  i3  an  S-3haped 
input/output  curve  which  is  typically  associated  with  hysteresis  and 
bistability.  A  rigorous  stability  analysis  shows  that  the  dotted  part 
(region  II)  may  correspond  to  an  unstable  steady  state,  resulting  in 
bistability  and  hysteresis,  as  shown  in  the  figure. 

6 .  Phase  Conjugation 

Optical  phase  conjugation  is  a  technique  that  employs  the  cubic 
nonlinearity  effect  of  the  optical  medium  to  precisely  reverse  both  the 
direction  and  the  overall  phase  factor  in  an  arbitrary  beam  of  light.  The 
process  may  be  thought  of  as  reflection  of  light  from  a  "mirror"  with 
unusual  image-transformation  properties.  In  order  to  illustrate  this  point, 
note  that  a  conventional  mirror  vfig.  24a)  changes  the  sign  of  the  k  vector 
component  normal  to  the  mirror  surface  while  leaving  the  tangential 
component  unchanged.  On  the  othr r  hand,  phase  conjugation  (fig.  24b)  cause 
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Fig.  23.  Hysteresis  curve  for  transmission  across  a  linear  nondispersive/ 
nonlinear  dispersive  interface  for  the  case  of  normal  incidence. 
The  nonlinearity  and  dispersion  parameters,  and  the  propagation 
constant  in  the  nonlinear  medium  have  exactly  the  same  values  as 
fig.  22;  and  the  ratio  of  the  (linear)  propagation  constants  of 
the  linear  and  nonlinear  media  is  taken  as  1/3. 


I 


reflection 
by  mirror 


"reflection' 
by  phase 
conjugation 


Fig.  24.  Principle  of  phase  conjugation  illustrated  through  a  comparison 
between  conventional  reflection  through  mirrors  and  reflection 
by  phase  conjugation. 
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an  inversion  of  the  entire  vector  k,  so  that  the  incident  ray  exactly 
returns  upon  itself  [39].  It  is  not  hard  to  argue  that  an  incident  converging 
(diverging)  beam  would  be  conjugated  into  a  diverging  (converging)  beam. 
Mathematically  speaking,  this  means  that  if  represents  the  incident 
phasor,  then  y*  would  represent  the  phase  conjugated  phasor. 


6 . 1  Comparison  with  Holography 

At  this  time  it  may  be  instructive  to  refresh  our  memory  with  the  basic 
principles  of  wavefront-reconstruction  imaging,  or  holography,  first 
proposed  by  Gabor  and  experimentally  realized  by  Leith  and  Upatneiks.  The 
wavefront-reconstruction  process  consists  of  two  basic  operations: 
recording  or  information  storage,  and  reconstruction  [19] .  To  understand  the 
first,  consider  two  wavefronts,  represented  by  their  respective  envelopes 


V.1,2  -  *1,2SXP  -j»l,2' 


to  interfere  in  space  (see  fig.  25a),  resulting  in  an  intensity 


1  '  '*.i  +  V 


lvall2+l\2|2*  *  \l\2 


2  2 

al  +  a2  +  2ai32  003  ' 


(140) 


The  intensity  distribution  is  then  usually  recorded  on  a  film  and  the  film 
developed  to  generate  a  transparency  or  "hologram"  with  a  transparency 
function  that  is  proportional  to  I.  To  understand  the  reconstruction 


(a)  (b) 


Fig.  25.  Imaging  by  wavefront  reconstruction  or  holography. 

(a)  recording,  (b)  reconstruction.  If  the  wave  illuminatin 
the  hologram  is  proportional  to  the  complex  conjugate  of  th 
reference  wave,  a  real  image  results. 
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process,  assume  that  thie  hologram  is  now  illuminated  with  a  reconstruction 
wave  Ve2'  amplitude  behind  the  transparency  is  then  proportional 

to 


WW2  *'V!*  VelVe2  * 


Note  that  if  y^  is  an  exact  duplication  of  one  of  the  original  waves,  viz. 

y  which  we  shall,  from  now  on,  call  the  reference  wave,  the  third  term  in 
e2 

the  above  expression  becomes  equal  to  y  multiplied  by  a  real  term,  and 
hence  is  a  duplicate  of  y  ,  which  we  shall  call  the  object  wave.  However, 
if  y^  i3  proportional  to  y^*,  observe  that  the  fourth  term  of  the  above 
expression,  which  we  shall  name  y^,  becomes  proportional  to  y^*,  and  hence 
represents  the  conjugate  of  the  object  wavefront.  We  may  state  at  this  point 
that  can  also  be  thought  of  as  being  responsible  for  the  formation  of  a 
real  image  of  the  original  object  (see  fig.  25b)  [19] . 

From  the  discussion  in  the  last  paragraph,  it  is  clear  that  one  can 
achieve  phase  conjugation  using  holographic  techniques,  although  in  the 
conventional  sense,  the  process  is  slow,  and  not  real-time,  owing  to  the 
efforts  in  making  the  hologram.  Note  that  in  the  above  discussion  we 
tacitly  talked  in  terms  of  envelopes  rather  than  phasors  since  reversal  of 
the  entire  k  vector  is  not  usually  possible  using  the  geometry  shown  in  fig. 
4.25.  The  process  of  introducing  y  ^  3  all  at  once  in  a  cubically 
nonlinear  medium  to  generate  y^4  speeds  up  the  conjugation,  and  is 
technically  referred  to  as  the  four-wave  mixing  problem.  The  entire  k 
vector  of  y^4  is  the  negative  of  that  of  y  .  This  will  be  discussed  at 


some  length  in  the  following  subsection. 
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6.2  Semiclasaical  Analysis 


We  start  from  eq.  (5)  with  P2  =  0: 


aVdt2  -  cVv  -  (2p3/3)82V3/3t2, 


(141) 


and  take  y  to  be  the  sum  of  four  waves  as 


V  “  Re[(eypl(z)  +  Vp2(Y/Z)  +  Vp3(y/Z)  +  eyp4(z))  exp  jcOQt]  (142) 


where  Vp2'  Vp3  represent  the  phasors  corresponding  to  the  (undepleted) 

"pump”  waves  with  y  ,  Vp4  representing  the  "probe"  and  the  "conjugate"  wave 

phasors  respectively.  The  quantity  £  is  taken  to  be  small  to  emphasize  that 

the  pump  energy  is  much  larger  than  the  probe  and  conjugate  wave  energies. 

Substitution  of  (142)  into  (141)  and  retaining  only  contributions  around  to 

o 

gives  the  following  equation  upon  equating  coefficients  of  E°: 


‘W  +  "o'W  '  <f3t^lVV2|VV' 


(143) 


where  kQ  =*  (iWc^.  Now,  assuming  the  pumps  to  be  contrapropagating  in  the 
sense  that  (see  fig.  26) 


V  “  *e  exp  j(kzz+kyy) 


V  ,  =  V  exp  +  j(k  z+k  y) ,  y  constant 
p3  e  z  y  e 


(144) 


and  of  equal  amplitude,  the  following  expression  for  the  "effective" 
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(NONDEPLETED  PUMP  WAVE) 


Fig.  26.  Basic  geometry  of  phase  conjugation  by  four-wave  mixing. 
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propagation  constant  k^  results  upon  substitution  in  (143) 


ko  ‘  (ky+k*,/t1‘(3 133/2)  ^e|2}- 


Now  equating  the  coefficients  of  £  yields 


(145) 


d2%l%4l/dz2  +  *'^p2.+'*,p4^ 


■|P3ko/2)  [l<lp2+Vp3|2<¥pl+VP4)  +  2IV+Vp3|2%1*VI 


(146) 


We  now  write 


Vpl  =  Vel  (z)exp  "  ^02' 


Vp4  =  Ve4(z)exp  +  jkQZ 


(147) 


in  anticipation  of  a  phase  conjugated  wave  y  traveling  in  the  -z  direction 

due  to  the  probe  wave  \y  (see  fig.  26) .  We  also  employ  the  slowly- 

varying  assumption  on  Vj/  and  \\i  .  Furthermore,  no  .*.•  that  the  term  (y  + 

ex  qh  pc 

2  2 

V  ‘"V  +  I  for  our  case)  defines  essentially  a  sinusoidal  phase 

grating  in  space  due  to  the  cubic  nonlinearity.  In  the  direction  of 

propagation  of  the  probe  and  the  conjugate  wave,  however,  only  the  term  from 

the  above  expression  that  is  spatially  constant  will  affect  their  evolution. 

The  spatial  distribution  of  the  induced  phase  grating  will  have  variations 

along  the  directions  of  propagation  of  the  pumps;  however,  these  will  not 

affect  the  probe  and  conjugate.  Mathematically,  this  amounts  to  picking  out 

2 

the  constant  term  from  (y^  +  V  •  EcI-  <146)  thus  gives  the  following  set 
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of  coupled  equations  between  y  and  y  : 


dyel/dz  -  2jp3kolVe!2<2yel+  ye*) 


dy  */dz  =  -2j(5,k  ly  l  (y  .  +  2y  *) 
e4  o  Te  Tei  e4 


To  solve  the  system  <148),  we  introduce 


Vei  =  Vei  6XP  "4jP3koIVe'  2 


y  -  y  .  exp  4 jp  k  ly  i  z 
e4  Te4  r3  o  e 


to  recast  the  system  in  the  form 


(148a) 


(148b) 


(149) 


dy  /dz 


dy  ./dz 
e4 


j8vei' 


(150) 


where 


8  -  2|i  k  |y  I2. 

3  o  e 


(151) 


The  solution  to  (150)  for  initial  conditions  y  , (0)  and  y  . (L)  is 

el  e4 

given  by  [40] 
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Vq1(z)  *  j  [  iSl  sin  (|5|z)/5  cos  ( |S  I  1>)  (L) 

+  (cos  ( |? | (z-L) ) /cos ( |S I L) ] V  (0) ,  (152) 

Ve4(z)  =  [cos  {  i§  I  z) /cos  (  |8  IL)  ]\j/e4  (L) 

-  j  [  I S  I  sin  (  |8  I  (z-L)  ) /§cos(  |8  IK)  ]\jjr  *  (0)  .  (153) 

In  the  practical  case  of  phase  conjugation,  \|/  (0)  is  finite,  whereas 

\y  (L) =  0 .  In  this  case,  the  nonlinearly  reflected  wave  at  the  input  plane 
e4 

(z=0)  is 

Ve4(0)  =  j  [  (8/  1 8  | )  tan  <|S|L)]v  J(0),  (154) 

where  8  is  defined  in  (151)  .  Note  that  Ve4  is  proportional  to  the  complex 
conjugate  of  as  expected. 

7 .  The  Nonlinear  Schrodinger  Ecruation  and  Soliton  Propagation 

It  is  well  known  in  the  theory  of  nonlinear  waves  that  pulses  traveling 
in  nonlinear  medium  distort  under  the  effect  of  nonlinearity  [1],  [2]. 

Furthermore,  from  linear  theory,  we  may  also  recall  that  pulses  distort  and 
spread  due  to  dispersion  in  the  medium  [4],  It  is  sometimes  possible,  in  a 
nonlinear  dispersive  environment,  to  ensure  distortionless  propagation  of 
pulses  as  a  result  of  a  balance  between  the  nonlinearity  and  dispersion.  For 
baseband  propagation,  it  is  a  little  easier  to  see  why  this  should  happen: 
the  steepening  effect  of  nonlinearity  (see  fig.  1)  may  be  balanced  by  the 
spreading,  or  "smoothening'"  effect  of  dispersion.  For  modulated  pulses,  both 
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nonlinearity  and  dispersion  cause  a  chirping  in  frequency  to  develop  during 
propagation,  as  well  as  other  distortions  of  the  envelope.  Given  the  right 
amounts  of  nonlinearity  and  dispersion,  it  is  possible,  onece  again,  to 
offset  the  effect  (s)  of  one  with  that  of  the  other  to  ensure  distortionless 
propagation  of  the  modulated  pulse  [2], 

In  keeping  with  the  nature  of  the  presentation  thus  far,  we  will  derive 
the  PDE  for  the  envelope  in  the  nonlinear  dispersive  medium  in  the  fastest 
possible  way.  Consider,  first,  a  linear  medium  with  arbitrary  dispersion, 
as  shown  in  fig.  27.  Since  we  are  discussing  the  behavior  of  modulated 
pulses,  it  is  probably  best  to  locate  ourselves  at  a  point  (co^k^)  of  the 
dispersion  curve  (the  coordinates  corresponding  to  the  carrier  frequency  and 
propagation  constant)  and  define  "excursions"  around  that  value  by  a  Taylor 
series  expansion  of  the  form  [2],  [4],  [41] 

co-co  =  u  (k-k  )  +  (u'/2)^,  co  >  0 .  (155) 

o  o  o  o 

2  2 

In  (155),  u  =  3(0/3k  |  and  u'  *  8  (0/3k  '  are  called  the  group  velocity 
o  co  o  co 

o  o 

and  group  velocity  dispersion  terms  respectively,  and  can  be  explicitly 
calculated  from  a  knowledge  of  the  explicit  equation  describing  the 
dispersion  curve  in  fig.  27.  Changing  to  variables 

=  co  -  co  ,  K  =  k  -  k  (156) 

o  o 

we  can  reformulate  the  nature  of  the  dispersion  around  the  carrrier 
frequency  C0q  as 

£1  -  u  K  +  (u1  /2)K2  . 
o  o 


(157) 
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How  do  we  figure  out  the  PDE  describing  envelope  propagation  ?  The 
procedure  to  be  followed  is  the  inverse  of  what  we  did  to  calculate  the 
dispersion  relation  (103)  from  the  linear  part  of  (101) .  Note  in  that  case, 
we  substituted  V  -  exp  j  (cot-kz)  in  the  linear  PDE,  which,  in  reality, 
amounted  to  replacing  3/3t  by  jco  and  3/3z  by  -jk,  where  0)  and  k  were  the 
respective  frequency  and  propagation  constant  parameters.  In  this  case,  we 
do  exactly  the  inverse:  we  replace  by  -j3/3t  and  K  by  jd/dz.  The 
wavefunction  these  operators  must  operate  on  is,  however,  the  envelope  y^  of 
y;  since  we  are  located  around  the  carrier  frequency  on  the  dispersion  curve. 
With  all  this,  (157)  yields  the  linear  PDE  [2] 

j(3\ye/3t  +  uQ3\)/e/3z)  -  (uV2)  32ye/3z2  =  0,  (157a) 

or 

j3y  /3t '  -  (u'/2)  d2y  /3z'2  =  0  (157b) 

e  o  e 

if  we  are  in  a  traveling  frame  of  reference 

z'  -  z-u  t;  t'  -  t.  (158) 

o 

In  what  follows,  we  will  modify  (157)  by  heuristically  incorporating 
the  effects  of  nonlinearity.  Assuming  a  cubic  nonlinearity  ((J^  *  0;  “  0) 

and  y  -  Re  {y  (z,t)  exp  j  (to  t  -k  z)}  in  (3b),  and  using  the  "slowly- 
varying  envelope"  approximation,  it  follows  that  yp  evolves  according  to 

2 

j3y  /3t'  +  (to  0  / 4 > \|/  y*  - 

©  O  J  ©  © 


o 


(159) 
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in  a  suitable  traveling  frame  of  reference.  In  our  heuristic  approach,  we 

will  include  the  effects  of  nonlinearity  in  (157)  by  simply  replacing  the 

2 

j3y  /3t'  in  (157b)  by  j3y  /3t'  +  ((o  P  /4,y  y*  .  This  yields 

©  6  O  O  ©  © 

j3y  /3t'  -  (u'/2)32y  /3z'^  +  (co  B,/4)y2y*  =  0,  (160) 

e  o  e  or3  e  e 

which  is  the  nonlinear  Schrodinger  (NLS)  equation  describing  envelope 
propagation  in  a  (cubically)  nonlinear,  dispersive  medium  [2],  [4],  [17], 

[20],  [41].  The  equation  has  been  extensively  used  to  model  waves  in  plasma 
[42],  water  waves  [43],  electron  and  ion  cyclotron  waves  [44]  and,  with  minor 
changes  of  the  independent  variables,  pulse  propagation  through  a  single  mode 
nonlinear  optical  fiber  [4] .  We  remark  that  the  nonlinearities  in  fibers 
may  be  Kerr-type,  while  the  dispersion  is  essentially  material  dispersion, 
since  no  modal  dispersion  effects  exist. 

To  find  particular  solutions  of  the  NLS  equation,  note  that  it  has  the 
same  form  as  eq  (71)  that  describes  beam  propagation  through  a  cubically 
nonlinear  medium.  In  fact,  an  interesting  analogy  of  pulse  propagation  in 
one  dimension  and  beam  propagation  may  be  made  by  comparing  eqns.  (160)  and 
(71) ;  in  fact,  diffraction  may  be  looked  upon  as  being  somewhat  like 
dispersion  in  space.  The  (real)  amplitude  of  a  particular  solution  is 
therefore  of  the  same  form  as  (75)  and  a  plot  appears  in  fig.  9,  with 
renamed  coordinates.  Let  us  conclude  by  saying  that  solitons,  as  these 
envelope  profiles  are  called,  have  been  genrerated  and  propagated  over  long- 
distances  through  optical  fibers  [21],  leaving  open  the  door  for  innovative 
communication  systems  that  can  be  much  faster  than  the  present  ones  that 
suffer  the  effects  of  material  dispersion.  The  advantages  of  solitons  are 
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that  they  not  only  remain  undistorted  by  themselves  in  a  nonlinear 
dispersive  medium,  but  can  also  interact  with  each  other  without  a  change  in 
shape  [20] . 

8 •  Conclusion 

This  report  reflects  some  of  the  work  done  or  being  done  in  the  area  of 
nonlinear  optics.  As  stated  in  the  Introduction,  it  is  by  no  means  complete 
or  exhaustive,  nevertheless,  I  hope  it  will  get  readers  started  in  the 
seemingly  complex  though  very  interesting  area.  Since  only  principles  have 
been  discussed  and  in  fair  generality,  the  discussion  could  find  some  use  in 
other  fields,  e.g.,  hydrodynamics  and  plasmas,  where  similar  effects  may  be 
observed  or,  or  least,  simulated.  References  are,  by  no  means  exhaustive, 
rather,  they  have  been  kept  at  a  minimum;  however,  interested  readers  can 
cross-reference  from  the  cited  work.  It  would  be  true  to  say  that  being  a 
starter  in  the  area,  I  could  appreciate  the  difficulties  one  could  incur, 
and  the  discussion  presented  above  is  just  a  summary  of  my  lines  of  thought . 
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